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Abstract: 

We present the ingredients necessary for the determination of physical K — > txtx decay 
amplitudes for AI = 3/2 transitions, from lattice simulations at unphysical kinematics 
and the use of chiral perturbation theory at next-to-leading order. In particular we derive 
the expressions for the matrix elements 7=2 ( im\Ow\K) , where Ow is one of the operators 
appearing in the AS* = 1 weak Hamiltonian, in terms of low-energy constants at next-to- 
leading order in the chiral expansion. The one-loop chiral corrections are evaluated for 
arbitrary masses and momenta, both in full QCD and in the quenched approximation. We 
also investigate the finite-volume effects in this procedure. 



PACS numbers: 11.15.Ha,12.38.Gc,12.15Ff 



1 Introduction 



The need for a quantitative control of strong-interaction effects in K — > tttt decays is 
underlined by the recent measurement of a non-zero value of e'/e [1] (demonstrating the 
existence of direct CP- violation) and the long-standing AI = 1/2 rule. A precise evaluation 
of non-perturbative QCD effects in nonleptonic kaon decays by lattice simulations, although 
possible in principle, is a truly extraordinary challenge, with a number of major theoretical 
and computational obstacles to be overcome. Ultimately, of course, we would wish to 
perform the simulations in full QCD, with physical quark masses on a lattice which is both 
large enough and sufficiently fine-grained for finite-volume effects and lattice artefacts 
to be negligible (or at least under control). We are some way from being able to do 
this, and therefore need rely on approximations. In this paper we explain our strategy, 
first presented in ref. [2], for improving one of the key approximations, the use of Chiral 
Perturbation Theory (%PT). Specifically we propose to perform the calculations at next- 
to-leading order (NLO) in the chiral expansion. In this paper we present the necessary 
formulae for AI = 3/2 decays and calculate the chiral logarithms for both full QCD and 
for the quenched theory. The chiral logarithms have been evaluated for arbitrary meson 
masses and momenta. 

In our calculations we incorporate the recent developments in our understanding of finite- 
volume effects in nonleptonic kaon decays [3, 4, 5]. At NLO in %PT, for AI = 3/2 
transitions, the corrections which vanish only as powers of the volume can be eliminated 
using the techniques presented in these papers, even in the quenched approximation (note 
however, that the generalization of the techniques in refs. [3, 4, 5] to the case when the two- 
pions have non-zero total momentum has yet to be performed). This point is explained in 
some detail in sec. 4 below. This is not the case however, when applying quenched %PT to 
AI = 1/2 transitions. For these decays, the lack of unitarity in the quenched theory leads 
to a number of subtleties and major difficulties which need to be overcome in order to 
understand the volume dependence and to extract the amplitudes [6, 7, 8]. We expect 
similar difficulties for AI = 3/2 transitions beyond one- loop order in chiral perturbation 
theory in quenched QCD. 

In this paper we do not discuss the ultra-violet problem, i.e. the construction of finite matrix 
elements of renormalized operators constructed from the bare lattice ones. For AI = 1/2 
transitions this involves the subtraction of terms which diverge as inverse powers of the 
lattice spacing (and for fermion actions which explicitly break chiral symmetry, such as the 
Wilson action, there is also mixing with other operators of dimension 6). The subtraction of 
power divergences is unavoidable, and different lattice formulations of QCD, including ones 
which preserve chiral symmetry, are being used in attempts to perform these subtractions 
effectively and to determine the matrix elements as precisely as possible. The strategy 
discussed in this paper for AI = 3/2 decays can be applied to any of these formulations of 
lattice fermions. 

The non-perturbative QCD effects in AI = 3/2 K — > tttt decays are contained in the 
matrix elements of the following three operators which appear in the AS = 1 effective 
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Hamiltonian: 



04 = {s a d a ) L (upup - d p dp) L + (s a u a ) L (updp) L ; (1) 
3 

Or = -(s a d a ) L e q(Wf3)R'i ( 2 ) 

q=u,d,s,c 

3 

08 = ^(s a dp) L J2 e q(ma)R, (3) 
q=u,d,s,c 

where a, f3 are colour indices and e q is the charge of quark q. The subscripts L, R signify 
"left" and "right" so that (qiq2)L,R = <?i7/u(l T 75)92- 04 transforms as the (27,1) repre- 
sentation of the SU(3)l x SU(3)r chiral symmetry and the electroweak penguin (EWP) 
operators 0-j$ transform as the (8,8) representation. O4 is a AI = 3/2 operator, whereas 
07,8 have both AI =1/2 and 3/2 components. The definition of the EWP operators 07, 8 
coincides with that used in refs. [9, 10] (where they were denoted by Qi^). Our definition 
of O4 is that introduced in ref. [11]. In terms of the operators defined in refs. [9, 10] O4 
is equal to the AI = 3/2 component of Qi + Q 2 . In this paper we study the K — > rnr 
matrix elements of these operators, with the two pions in an / = 2 state, so that only the 
AI = 3/2 components of (9 7;8 contribute. 

In this project we are evaluating K — > nn matrix elements directly. All but one of the 
previous lattice studies evaluated K — > tx matrix elements of the form ( tt\Ow\K) and used 
lowest order to relate them to matrix elements ( tttt\Ow\K), following the original 

proposal in ref. [12] (see also refs. [13, 14]). (CV represents any of the AS = 1 operators 
which appear in the effective weak Hamiltonian.) The exception is the quenched study of 
the CP-conserving K — > nn AI = 3/2 decays in ref. [15], where the two pions were at rest. 
The previous studies are therefore not sensitive to final-state interactions. 

In ref. [2] it was shown that it is possible to obtain all the low-energy constants necessary 
to determine the physical A J = 3/2 decay amplitudes at NLO in the chiral expansion 
by performing simulations with a simple range of momenta for the mesons 1 . Specifically 
the kaon and one of the final-state pions are at rest and the second pion has an arbitrary 
momentum (the energy of this second pion will be referred to as E n in the following). In 
order to ensure that the final state has isospin equal to 2, we have to symmetrise the final 
state over the two pions. Since K — > tttt matrix elements with this set of momenta are 
being computed by the SPQ C( jR collaboration [16], we refer to it as SPQR kinematics. Of 
course, it is also possible to determine the low-energy constants which appear at NLO in 
the chiral expansion from simulations using other kinematical configurations 2 . For this 
reason in the following we will present the results with arbitrary kaon and pion masses and 
momenta. 

Among the previous lattice studies are two recent determinations of K — > mr matrix ele- 
ments from computations of K — > n matrix elements with quenched domain-wall fermions 

lr The generalization of the formalism to AI =1/2 decays is currently under investigation 
2 In ref. [17] it is suggested that an alternative set of kaon matrix elements be used to determine the 
low energy constants. We comment on this proposal in sec. 6. 
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and using %PT at lowest order [18, 19]. The very interesting results for the octet enhance- 
ment, ReA /ReA 2 = 9-12 [18] and 24-27 [19] (compared to the experimental value of 
22.2) and for e'/e= {(-7) - (-2)}1(T 4 [18] and {(-8) - (-4)}1(T 4 [19] (compared to the 
experimental value of +(17.2 ± 1.8) 1CT 4 ), will serve as important benchmarks for future 
calculations. (Ai is the K — > nit amplitude for decays into two-pion states with isospin I.) 
These papers also contain detailed discussions of the systematic errors and the difficulties 
in estimating and controlling them. The importance of a reliable error analysis is under- 
lined by a comparison of the results of the two collaborations and with the experimental 
values. For example, the cancellations which occur between different contributions to e'/e 
(in particular between the matrix elements of operators conventionally called Oq and 0$) 
imply that each component must be determined with good precision. In this paper we 
tackle one of the main sources of systematic error, the use of %PT at lowest order. 

The remainder of the paper is organized as follows. In the next section we explain our 
strategy in detail. We express the matrix elements in terms of the low-energy constants 
and demonstrate that it is possible to determine the physical matrix elements at NLO in 
the chiral expansion from numerical simulations which are currently feasible. We present 
a brief introduction to quenched and unquenched xPT in sec. 3. Sec. 4 is devoted to a 
study of finite volume effects in %PT . The NLO results are presented in sec. 5 for physical 
kinematics, i.e. with the momenta of the mesons satisfying p K + = p n + + p n o, but with 
the masses being arbitrary. In this section we also present the results for the special cases 
itlk — ^va^ and itlk = rn n , in both cases with each of the pions at rest, and compare 
our results with those of previous papers where this is possible. The expressions for the 
matrix elements at NLO in the chiral expansion for the SPQR kinematics are presented in 
appendices B.l and C.l for full QCD and the quenched approximation respectively. The 
results for general kinematics, i.e. with arbitrary masses and momenta, are very lengthy 
and we present them on the web site [23], from which interested users can import them 
into their programs. For full QCD we also present the results for general kinematics in 
appendix B.2. Sec. 6 contains our conclusions. There is one other appendix which contains 
integrals used in the study of finite-volume effects. 



2 The Strategy 

In this section we present our strategy for the determination of K — > nn matrix elements 
using chiral perturbation theory at next-to-leading order. We envisage computing the 
matrix elements in numerical simulations for a variety of (unphysical) masses and momenta. 
The generic form for the behaviour of the matrix element of a weak operator Ow with 
masses and momenta is 

( tt+tt \Ow | K + ) = Ap(l + chiral logs) + ]T A^, (4) 

i 

where A and the Bi are kinematic factors and p and the Aj are low-energy constants. (In 
the applications below p corresponds to c^ 27 ' 1 ) or 7( 8 ' 8 ) and the Aj's correspond to the /3j's 
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or the 5^.) A is 0(p 2 ) when O = 4 and 0(p°) when O = C 7 ,s and the B t are 0(p 4 ) 
and 0(p 2 ) when (9 = 4 and (9 = 7:S respectively. In both cases there is only one low- 
energy constant at lowest order of x?T and several at next-to-leading order. For the EWP 
operators the low energy constants p and A, include a factor of the electromagnetic 
coupling e 2 . The expression "chiral logs" on the right-hand side of eq. (4) is shorthand for 
the full one-loop contribution, which includes both terms containing logarithms and also 
others without logarithms. 

Of course we wish to determine the matrix element for a physical kaon decaying into two 
physical pions. Our strategy for doing this at NLO in the chiral expansion is as follows [2]: 

1. Calculate the chiral logarithms at one-loop order in for general masses and 
momenta. In this paper we perform this calculation for the AI = 3/2 operators 
04,7,8 in both full QCD and in the quenched approximation. Once the logarithms 
have been computed, we have the formula for the behaviour of each matrix element 
in terms of the low-energy constants p and Aj. 

2. Determine the matrix elements in eq. (4) using lattice simulations for a variety of 
masses and momenta. The kinematic range must be sufficient to enable the low- 
energy constants to be obtained from a fit to the behaviour of the matrix elements 
with the masses and momenta. 

3. Use the low-energy constants determined in item 2 and the expression in eq. (4) with 
itlk and m n equal to their physical values and zero energy-momentum transfer at 
the operator (i.e. with p K + = p n + + p n o) to determine the physical matrix elements. 
In extracting the physical amplitudes we assume that the higher order terms in the 
chiral expansion can be safely neglected in the region of the extrapolation. 

For compactness of notation, in the following we represent the low-energy constants by the 
same symbols in both full QCD and in the quenched theory. It must be stressed however, 
that the low-energy constants are not expected to be the same in the two cases. Indeed 
their dependence on the renormalization scale is in general different. 

In the remainder of this section we give the expressions for the matrix elements of 04,7,8 
in terms of the low energy constants. 



2.1 Counterterms for O4 up to 0(p 4 ) 

The general expression for the matrix element of O4, which transforms as a (27,1) under 
the SU (3) x SU (3) chiral symmetry, for arbitrary meson masses and momenta is 



(7r + 7r |C 4 |^ + ) = 




cr ' ( 8^0 • p K + 6p w o ■ p n + — 2p n + ■ p K + chiral logs 



+(3 2 ( -24mlm 2 K + 24ml) 
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+/3 4 ( 16m 2 (p n o -p K ) + Vom\ (p n o ■ p K ) + 24m 2 (Ptt° • Pw+) 
-4m 2 (Pn+ -Pk)- 4m^ (p w + -p K )) 
+[3 5 (-16m 2 (Ptt° • p x ) - 12m^ (p^o • p u +) + 4ml (Pn+ • Pa-)) 
+ / 9 7 (l6m^ (p^o • p x ) + 32m^ (p^o • p x ) + 12m 2 (p n o ■ p 7r +) 

+24m 2 K (p^o • p n +) - Ami (Pn+ • Px) - 8m^ (p 7r + • p^)) 

+/3 22 (24(p X • Pno) {PK • Ptt+) - 8(PA- • p^o) (p^o • p^) + 32 (j) K ■ p w +) (p n •£>„■+)) 

+/3 24 (l28(p i , -p^o) 2 - 96(p^o -p^ + ) 2 - 32(p x -p^ + ) 2 ) } • 



(5) 



/tt and fx are the physical pion and kaon decay constants. At NLO in the chiral expansion 
they are given in terms of the parameter / of the chiral Lagrangian (see sec. 3 below) and 
the Gasser-Leutwyler coefficients in ref. [20] for full QCD and in ref. [21] for quenched QCD. 
a (27,i) j s tu e leading order (0(p 2 )) low-energy constant whereas the $ (i = 2, 4, 5, 7, 22, 24) 
are the 0(p 4 ) counterterms, where the subscript corresponds to the numbering of the 
operators in ref. [22]. These operators are listed explicitly in eq. (13) of sec. 3. Thus 
for a set of simulations, with a range of quark masses and momenta in the chiral regime, 
the mass and momentum dependence is given by eq. (5). We envisage determining the 
low-energy constants a^ 27,1 ) and $ in lattice simulations (or at least a sufficient subset of 
these constants). We have calculated the chiral logarithms for arbitrary meson masses and 
momenta, for full QCD and in the quenched approximation and we present the results on 
the web site [23]. For full QCD we also present the results in appendix B.2. 

It is possible to determine the physical matrix elements of O4 by performing simulations 
with the SPQR kinematics. This will become clear in the following subsections where 
we rewrite eq. (5) for physical kinematics (subsection 2.1.1) and for SPQR kinematics 
(subsection 2.1.2). The corresponding one-loop expressions with the chiral logarithms 
included are given in: 

• eqs. (79) - (82) for the physical kinematics in full QCD; 

• eqs. (91) - (94) for the physical kinematics in quenched QCD; 

• eqs. (108) - (116) for the SPQR kinematics in full QCD; 

• eqs. (133) - (136) for the SPQR kinematics in quenched QCD. 

2.1.1 The Matrix Element of 4 for Physical Kinematics 

For the physical matrix element, when there is no momentum insertion at the weak operator 
04, eq. (5) reduces to 
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+(ft - ft + 4ft + 2(3 22 )m% + (4ft - 4ft - 2ft - 16/3 24 )m^ 

+(-4ft + 3ft + ft - 2ft ~ 2ft 2 + l6(3 2 4)m 2 K ml j . (6) 

Of course it must be remembered that the chiral logarithms depend on the kinematics, 
however, as they are calculable we do not exhibit this dependence explicitly here. In 
principle, by varying the masses and momenta of the pions it is possible to determine 
c^ 27,1 ) and the ft's but in practice this is not the optimal strategy For each set of pion 
masses and momenta one would have to tune the mass of the strange quark in order to 
ensure energy conservation, and in addition one would have to extract the matrix elements 
of non-leading (excited) energy levels. In the following subsection we demonstrate that 
by performing the simulations with the SPQR kinematics it is possible to determine the 
combinations of the low-energy constants which appear in eq. (6). 

2.1.2 The Matrix Element of 4 for the SPQR Kinematics 

For the SPQR kinematics the matrix element in eq. (5) reduces to 

(7r + 7r |C 4 |^ + )sPQR = -p^la (27 ' 1) ( K E n m n + ^m K (E 7T + m w ) + chiral logs) 

+4ftr< + (4ft + 2fr)E n ml + (ft - ft + fc)m\m K 
+ (ft - ft + ft + 2(3 22 )E w mlm K + (-4ft + 8ft 4 )m> 2 : + (ft + 2(3 7 )E n m 3 K 
+(-2/3 5 + 4ft + 4(3 22 )E n m n m 2 K + (ft + 2(3 7 )m n m z K 

+ (-16ft 4 )£ 2 m 2 + 2ft 2 £> 7r m K + 8(3 24 E 2 w m 2 K j . (7) 

It can readily be seen that by fitting the lattice results for the matrix elements as a function 
of the masses and E n , the combinations of low-energy constants needed to determine the 
physical matrix elements in eq. (6) can be obtained. 

2.2 Counterterms for Oj j8 up to 0(p 2 ) 

We now present the corresponding formulae for the AI = 3/2 EWP operators 7>8 . In 
this case the matrix elements start at O(p ) and at NLO we have to consider all possible 
terms of 0{p 2 ). The general expression for these matrix elements up to this order is 

(n + 7T \O 7>8 \K + ) = | 7 (8 ' 8) (1 + chiral logs ) + 

JK J w { 

<*1 (3 (Pit ■ Vk) - (Ptt° ■ Ptt+) + g (Pir+ ' Ptf)) ~ 5 2 ((j^O • p K ) + (p w ■ p w + ) ) (8) 

S3 (^(p w o • p K ) + g(p ff + • P*)) + (^4 + 5 5 ) (2m 2 , + 4m 2 ) + 5 6 (4m 2 , + 2m 2 ) 1 . 
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The 0{p 2 ) counterterms 5i, i — 1 - 6, correspond to the operators with the same numbers 
in ref. [24]. We give these operators explicitly in eq. (21) of sec. 3. A complete list of the 
0{p 2 ) independent counterterms can be found in ref. [25]. 

Since O-j and 0$ both transform under the same representation, the generic form of the 
chiral expansion is the same but the values of the 7( 8 ' 8 ) and the Si on the right-hand side 
of eq. (8) depends on the operator. We leave this dependence implicit. 

We present the results for the chiral logarithms for arbitrary meson masses and momenta, 
together with those for C 4 , on the web site [23]. For full QCD we also present the results 
in appendix B.2. 

We rewrite eq. (5) for physical and SPQR kinematics in sec. 2.2.1 and 2.2.2 respectively. 
The one-loop expressions with the chiral logarithms included are given in: 

• eqs. (87) and (88) for the physical kinematics in full QCD; 

• eqs. (95) and (96) for the physical kinematics in quenched QCD; 

• eqs. (117) - (120) for the SPQR kinematics in full QCD; 

• eqs. (137) - (140) for the SPQR kinematics in quenched QCD. 

2.2.1 The Matrix Element of 7 8 for Physical Kinematics 

For the physical kinematics eq. (8) reduces to 

(n+n \O 7 , 8 \K + ) phjs = | 7 (8 ' 8) (1 + chiral logs ) + (9) 

( - 5 2 - 5 3 + 2(5 4 + 5 5 ) + 4<5 6 ) m\ + (b x + 5 2 + 4(5 4 + 5 5 ) + 2<5 6 )m' j. 

2.2.2 The Matrix Element of O t , 8 for the SPQR Kinematics 

Following the same strategy as for 4 , we envisage determining the 7( 8,8 )'s and <5;'s by 
performing the simulations with the SPQR kinematics. In this case eq. (8) reduces to 

( tt+tt \0 7>8 I K+ ) SPQR = | 7 (8 ' 8) (1 + chiral logs ) + 2(6 4 + 5 5 + 25 6 ) m\ + 

JK U { V 7 

_ s 2 - 25 3 )(m 7r + E v )m K + 2 (2(5 4 + 8 5 ) + 5 6 )ml - + 5 2 )m^ j. (10) 

It is straightforward to verify that if one determines 7( 8 ' 8 ) and the combinations of the <5j's 
which appear in eq. (10) by fitting lattice data for SPQR kinematics, then one has all the 
necessary ingredients to determine the physical matrix elements using eq. (9). 
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3 Quenched and Unquenched %PT for K —> tttt Decays 



In this section we outline the main ingredients of the one-loop xPT calculations. The 
generic form of the (27,1) component of the weak Lagrangian for K — > nn up to 0(p 4 ) 
is [22] 

C™ = -a™0™- E A0£' 1) +h.c. + O(p 8 ) (11) 

i=2,4,5,7,22,24 



where O^ 7 ' 1 ^ is the leading order (0(p 2 ) ) operator 



and O^J' 1 ^ are the operators which appear at NLO (0(p 4 )) 

OfV ] = -7S'(P)* < (P) , ,-, 



'4,2 

^4,4 — 

/ot 27 ' 1 ) 
^4,5 

/n(27,l) _ 

^4,7 _ 

/n(27,l) _ 

^4,22 — 

z^ 27 ' 1 ) - 

^4,24 — 



ruLM^, sy j: 



The tensor is defined as 



nl3 
12 



r 23 
■22 



Tn31 

1 12 



i32 



22 



Th13 

J 21 " 
1 

~2' 



Th31 

J 21 



1 

2' 



with all other components zero. L M and are defined by 

L M = i S f <9 M E, 



and £ and P by 
with 



■ iiv — 2 [d^Ly + OyL^) , 

5 = st x + x t s and P = ^( S t x - x t S ) 
-2(0|uu + dd|0) 



/ 5 



(12) 



(13) 



(14) 



(15) 
(16) 
(17) 



£ is the conventional exponential representation of the pseudo-Goldstone boson field 

r + K + \ 



cxp 



2i$ 

7T. 



( + JL 



with $ 



7T 



7T 



V2 + Ve n 



K° 



(18) 



f 7 ? J 



and M. is the quark mass-matrix M. = diag(m u , ma, m s ). 

For the (8,8) component of the AS = 1 weak Lagrangian the generic form is [24, 25]: 

£g) = - 7 (w)or-E^og 8 ) + o(p 4 ), 

1=1 

where O^ 8 ' 8 ^ is the leading order (0(p°) ) operator 



C^' e; = tr 



A 6 St Q E 



and are the NLO (0(p 2 )) operators 



^2,1 



o 

o. 
o 
o 



2.2 

(8,8 
2,3 



2.1 



2.5 



2,6 



= tr [A 6 

= tr [A 6 Lfj] tr [s f QSL M 
= tr [A 6 {StgS,L^Ln 
= tr [Ae^QE,^}' 
A 6 [s f Q E, P 
AgEtQE] tr[S] 



i tr 
tr 



Q and A6 are the charge and Gell-Mann matrices respectively: 



Q 




and A6 



/ 
1 
V 1 



(19) 



(20) 



(21) 



(22) 



For the strong interactions we only need the leading order chiral Lagrangian 

P 



■^-strong - tr 



(23) 



The NLO terms in £ s trong, which are proportional to the Gasser-Leutwyler Li coefficients, 
contribute to the wave function renormalization. These effects are reabsorbed in the renor- 
malization of the decay constant leading to the replacement of the factor 1/ f 3 by l/(f K p) 
in all the K — > nn matrix elements. We therefore do not present the NLO terms in the 
strong chiral Lagrangian. 



3.1 Quenched xPT 

For the calculation in quenched we use the strong Lagrangian introduced in ref. [21] 
P 



•-strong 
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-str 



(9 M Et)(^E)+Et x + x t S ] -ml^ + a (d,%)(d^ ) . (24) 
9 



In eq. (24) the trace over the chiral group indices in eq. (23) has been replaced by the 
supertrace over the indices of the graded group SU(3\3) L x SU(3\3) R . The fields £ and x 
are now graded extensions of £ and \ defined in eqs. (18) and (17). Since in the quenched 
approximation the singlet field remains degenerate with the pion and has to be treated 
as dynamical, the graded extension of the matrix $ is not supertraceless. We define the 
singlet field $ =str[$]/v / 6 proportional to the graded extension of the rj'. 

In the weak operators each field in eqs. (12), (13), (20) and (21) is replaced by its graded 
extension, and every trace by the corresponding supertrace. Although we use the following 
extension for the charge matrix in the quenched approximation 



the results for AI = 3/2 transitions do not depend on the choice of extension in the SU(2) 
(isospin) limit. This is because the EWP operators are always inserted on external quark 
lines, and hence the matrix elements are independent of the charges assigned to the ghost 
particles. For the same reason, for the present calculation we can simply extend A 6 to a 
six-by-six block diagonal matrix: 



In this section we study, at one-loop order in , the Euclidean correlation functions 
relevant for the finite-volume lattice determination of K — > nn decay amplitudes. As 
argued in the section 2, such a study is a necessary step to extract, by extrapolation, the 
physical amplitudes from correlation functions which at present can only be computed at 
unphysical values of masses and momenta [2]. Although the explicit discussion is presented 
for full QCD, it also applies to quenched (and partially quenched) QCD for AI = 3/2 decays 
(for AI =1/2 decays in quenched QCD this is not the case [8]). 

A significant step towards enabling the determination of the physical matrix element of 
(9 4 from matrix elements computed in lattice simulations at unphysical kinematics, has 
been the evaluation, at one-loop order in perturbation theory in finite volumes, of K — > nn 
amplitudes with mx = rn n and uik = 2m n and with all the mesons at rest [26, 27, 28]. 
(For a discussion of the advantages of using these kinematic points in the evaluation of 
AI = 1/2 decay amplitudes at leading order in the chiral expansion see refs. [29, 30].) The 
one-loop generating functional in the presence of weak interactions, also in the quenched 
approximation, was considered in ref. [31]. 

In discussing the behaviour of the correlation functions and matrix elements with the 
volume, we envisage varying the volume at fixed physics, i.e. with the masses, energies 
and momenta constant in physical units. From our explicit computations we see that all 
the one-loop corrections can be classified into five categories: 



Q = diag(2/3, -1/3, -1/3, 2/3, -1/3, -1/3), 



(25) 




(26) 



4 %PT on a Finite Volume 
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1. The infinite- volume chiral corrections to the matrix elements of the weak operators. 
These diverge in xPT , and are renormalized by operators of higher dimension [22] , 
see sec. 5; 

2. The finite- volume shift to the two-pion energy, AW = W — 2E W , which is independent 
of the weak operators, and only depends on the isospin of the two-pion state. We 
illustrate this point by showing explicitly that our calculations give the same shift 

3/2 3 /2 

AW for the operators O4, O r and O s (where the superscript 3/2 indicates the 
AI = 3/2 component of the operators); 

3. The corrections which shift the energy argument (E) of the amplitude A(E) from 
its infinite volume value to the finite volume one, i.e. A(2E 7T ) — > ^4(W). We show 
that in earlier studies the shift of the argument of A(E) was wrongly interpreted as 
a finite- volume correction to the amplitude [26, 27]; 

4. The finite volume corrections to the sink (source) which is used to annihilate (create) 
the two pions. These corrections obviously depend on the sink/source and can be 
eliminated non-perturbatively by studying the source-sink correlation functions. In 
the following we consider two classes of correlation functions. In the first the two 
pions are annihilated at the same time (as used, for example, in ref. [3]) and in the 
second one pion is annihilated much later than the other (as used, for example, in 
refs. [4, 5, 32]). 

5. Finally we have the Lellouch-Liischer (LL) finite- volume corrections to the matrix 
elements themselves. 

A very important result is that the finite-volume correction to the amplitude, in one-loop 
%PT and in the center of mass frame, is the universal LL factor with both types of sources 
mentioned in item 4 above. This will be explained in sec. 4.1 below. The generalization 
of the LL factor to the case where the two pions have a non-zero total three-momentum is 
being studied (see sec. 4.1.2). 

The starting point for our studies is the Euclidean correlation function 

C^(t K ,h,t 2 ) = (01^(^)^(^)0^(0)4(^)10) , (27) 

where 

Tr^t) = J d 3 x ir(x , t) e*® (28) 

is the Fourier transform of the one-pion interpolating field (with a similar definition for 
the kaon) and tx < < t% < t±. Ow is one of the AI = 3/2 weak operators, and from 
the correlation functions in eq. (27) we wish to determine the matrix element {nn\Ow\K). 
In the following subsections we explain, at one-loop order in xPT , what the finite-volume 
effects in the correlation function are, and how one recovers the LL factor in the matrix 
elements. 
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In order to simplify the discussion of the different corrections occurring at one loop, while 
maintaining all the essential ingredients, we begin in sec. 4.1 with a study of the correla- 
tion functions in a A$ 4 theory (for the strong interactions) and later, in sec. 4.2 and 4.3 
generalise the discussion to x?T . 



4.1 K + — > 7r + 7r Decays on a Finite Volume in a A$ Model 



In this section we illustrate the structure of the relations between correlation functions in 
finite and infinite volumes, by considering K — > ixix decays in a simple model in which the 
interactions between the pions are described by a A$ 4 theory. The bare K — > ixix vertex is 
taken to be a constant {\k-k-k)- This model contains all the essential ingredients necessary 
to illustrate the key features of finite- volume effects. In the following section we generalise 
this discussion to %PT . 

At tree-level, the correlation function in eq. (27) is given by 



CfUt K MM) = A^\E T ) x Q ( -°\t K ,t 1 ,t2,E T ) 



Xktttt X 



1 



2m K 2E& 2E& 



. e -m K \t K \ g-B^ti-ta) e -E T t 2 



(29) 



J Q2 



where E% = 

HI 



q-i\ + m 2 n and the total two-pion energy is given by E T = E^ + E$ 2 . 
A^(E T ) (= Xktttt) is the tree-level component of the K — > rrn decay amplitude. 

At one-loop order, we insert an interaction vertex (O s = —A |7T°(a; s )| 2 \n + (x s )\ 2 /A, where 
Ti represents the pion field) into the correlation function, and its contribution is given by 
the expression 



r[7r fl (ti)7ra(* 2 )CW0) 



/ 



d x s O s (x s ) 



Kl(t K )]\0) (30) 



corresponding to the Feynman diagram in fig. 1(b). In this case the weak operator is 
simply given by the triple scalar vertex, O w (0) = \ K7T7T K(0) ^(0) 7r^(0). 

In x?T , of the four diagrams shown in fig. 1 (and the graphs which renormalize the 
external meson fields) it is only diagram (b) which contains finite-volume effects which 
decrease as inverse powers of the volume. For the remaining diagrams the finite volume 
effects decrease exponentially and will be neglected. It is for this reason that the evaluation 
of this graph in the A<3> 4 model illustrates the key properties of finite volume effects in x?T . 
Note however, that the evaluation of the graph in fig. 1(b) is different in this model and 
%PT (for example both the weak and strong vertices in %PT depend in general on the 
momenta). 

The contribution from the diagram of fig. 1 (b) to the correlation function is: 



Xktttt A 



e -m K \t K \ Y 

2m K 1? 



I 



d,t 



g— 1*1— ts\ e —Etf 2 \t2—t s \ 



2E, 



91 



IE, 



k,l 



dkd did 



-i{k4,+U)t a 



-S 3 - - 

2tt 2tt {kl + wl){ll + wl) k+l ^ 
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Figure 1: One-Loop Chiral Perturbation Theory Diagrams for the K + — > 7r°7r + decay. The 
grey circle (square) represents weak (strong) vertices. In addition one must also calculate 
the contributions from the diagrams which renormalize the wave functions of the external 
mesons. 



— ^Ktttt A x 



e -m K \t K \ ^ r g _ -%il*l— *s| e ~Eq 2 \t2-t s \ g— U!l|ta| e -m\ts\ 



2m K L 3 ^ k + l >QJ 2Eff, 2E* 2 2 Wl 2w 2 

k.l 



e -m K \t K \ II I 

= -A fc Ax- — EfnV {T 1 + T 2 + T 3 + T 4 }, (31) 

where 

w \ = mj + \k\ 2 , w 2 = m 2 + | Z "J 2 (3 2 ) 

and <f = <fi + g*2- Xi— T4 are the contributions to the integrations over t s from the intervals 
(—oo,0), (0, t 2 ), (*2)*i) an d (ti,oo) respectively and are given by 



Eft + Eft +w 1 + w 2 

Eft + Eft - (Wt + W 2 ) 



p -E q - 2 (t 1 -t 2 )-(w 1 +w 2 )t 1 _ -E ?1 (t 1 -t 2 )~(w 1 +w 2 )t 2 

Ts = = =— ^ (33) 

Eft - Eft -w x -w 2 
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T 4 



e -Ej 2 (t 1 -t 2 )-(wi+w 2 )t 1 

Eft + Eft + w 1 + w 2 



For simplicity we choose \q\\ = \q 2 \ and take m\ = m 2 , so that q = 0, W\ = w 2 = w 
and E 1 ^ = Eft = E. The generalization to the SPQR kinematics is straightforward in 
perturbation theory, but as described in sec. 4.1.2, the generalization of the LL-factor to 
the case where the two pions have a non-zero total momentum has not been performed 
yet. Factorizing the tree level correlation function, we obtain 

C { }l{t K ,t u t 2 ) = -g^(t K ,t 1 ,t 2 ,2E) x A x (34) 

1 A ( -°\2w) f 1 e 2(E-w)t 2 _ I e 2Et 2 -2wt x _ e 2(E-w)t 2 e 2Et 2 -2wt 1 ~\ 

[2w) 2 \2(E + w) + 2{E-w) + ^2w + 2[E + w) J ' 

Note that in this specific example the weak amplitude is independent of the energy and 
thus A^(2w) = \ktt-k- The four terms in braces in eq. (34) are 7\ - T 4 respectively. For 
illustration we now consider the two cases < t\ = t 2 = t [3] and < t 2 t\ [4, 5, 32]. In 
the first case, T 3 = and we have a factor 

1 1 I e 2 ^-^ 1 e 2 ^-^ 1 - 1 1 \ 

1? \ J2wf \2(E + w) + 2{E-w) + 2{E + w)) ' (35) 

In the second case, since t\ ^> t 2 we can neglect T 4 since it is exponentially suppressed and 
the corresponding factor is 

1 v 1 j e 2 ^ E - w ^ e 2( - E - w ^ - 1 1 \ 

IJ^J^X 2w + 2{E-w) + 2(E + w))' ( ' 

Here we are considering finite volume corrections to correlation functions. As anticipated 
above, we will show that the finite volume corrections to the amplitudes, extracted after 
dividing by the appropriate source factor, are identical in the two cases, in spite of the 
differences of the expressions in eqs. (35) and (36). 

— * 

We now consider the sum over the momentum k and separate the contribution correspond- 
ing to w = E, with the two pions in the loop on the mass-shell, from the remaining terms. 
For the two cases we are studying, the contribution with w = E is given by 

— ^^/— +t+±=) (t 1 = t 2 =t) (37) 



and 



L 3 (2E) 2 VAE AE) 
f 1 1 1 

J 1- f n -I \ 



LH2E)n2E +k+ 4ES (U>>t ^' (38) 

where v = J2k- W= E- The term linear in time in eqs. (37) and (38), after resummation, 
corresponds to a shift of the two-pion energy due to the A<3> 4 interactions: 

W = E T + AW = 2E + . (39) 
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It is straightforward to check that the energy shift in eq. (39) corresponds to that obtained 
from the 0(A) term in the perturbative expansion of the Liischer quantization condition [33, 
34] using the one-loop expression for the s-wave phase-shift. 

Of the remaining terms in each of eqs. (37) and (38), the first one comes from final-state 
pion rescattering and will be cancelled by dividing by the matrix element of the source 
(this will be demonstrated explicitly below). These final-state interactions are different in 
the two cases ti = t 2 and t 2 <C t±. The third term in each of eqs. (37) and (38) is part of 
the (finite volume) one-loop correction to the "weak" matrix element and, as expected, is 
the same in both cases, t 1 = t 2 and t 2 <C t\. 

It is convenient to group the terms with w ^ E as follows: 

h ^TkAj^f^' + ^} (tl = h = t) (4o) 

In the limit t 2 — > oo, the terms proportional to e 2 ^ E ~ w ^ t2 with w > E are exponentially 
suppressed and only a finite number of contributions with w < E survive. As pointed out 
by Maiani and Testa, these terms dominate the correlation function over the contribution 
one is trying to isolate, i.e. the one with w = E. We denote these terms as Maiani-Testa 
terms. They can be eliminated by using the behaviour of the correlator as a function of 
t 2 (at large t 2 ) and separating the different exponentials, corresponding to the different 
energies [3, 4, 35]. We assume that this can be done in the numerical simulations and do 
not consider these contributions anymore. Up to exponentially small terms, the final term 
in eqs. (40) and (41) can be written as [34] 

— V 1 W = f( F ) | Zn | V (AO) 

L 3 ~ (2w) 2 w 2 -E 2 J{ ' 4(2n) 2 EL 8E 3 L 3 ' 1 } 



where n is defined by <fi = ^n, 



— I f°° d 3 h 1 ] 

= E p2 ^ . and m = V 7^ r ^' (43) 
,-, (i — n z ) J-oo 2tt Y Aw w z — E l 



\A £ \n\ 



(2tt) 3 Aw w 2 - E 2 



It is instructive to rewrite f(E) in a form which makes clear the connection to the one-loop 
contribution to the physical K — > -kit amplitude. We write f(E) =Re{—ifM(E)} where 

f M (E) = J j—^ ^ 2 _ m2 + i£ ^ q _ k y_ m2 + i£ ^ ( 44 ) 

is the relevant integral in Minkowski space (hence the subscript M). The relation f(E) = 
Re{—ifM(E)} can be readily verified using contour integration in /c -space. The real part 
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of the K — > 7T7T amplitude in Minkowski space is proportional to Re{l + iXfui^E)} = 
1 — Xf(E), so that one correctly obtains the real part of the Minkowski amplitude from 
the Euclidean correlation function, in agreement with the general arguments of ref. [32]. 
(At this order in perturbation theory the modulus of the amplitude is equal to its real 
part.) We can also rewrite f(E) as a one-loop four dimensional integral in Euclidean 
space, f (E)=Re{f E {iE)} , where 



Ref E (iE) = Re 



d A k 



1 



(2tt) 4 (k 2 + m 2 )((q - k) 2 + m 




(45) 



-AE 2 



After the integration has been performed in Euclidean space, the substitution q 2 = —AE 2 
is made. 

From the above discussion we conclude that the correlation function, at one-loop order in 
perturbation theory and keeping only the terms with the required dependence on t 2 , can 
be written in the form (for the two cases being considered): 



V*i(**,M) = A {0 \W)xg^(t K ,t,t,W)x 

3 



1 - A 



Ref E (iW) 



V 



4(2tt) 2 EL L 3 1QE 3 L 3 1QE 3 



and 



c^{t K MM = AW(w)xgM(t K ,t 1 ,t 2 ,w)x 



A 



Ref E (iW) + 



+ 



v 



+ 



A{2n) 2 EL L 3 1QE 3 L 3 8E 3 



(46) 

-0(\ 2 ) (t 1 = t 2 = t), 
(47) 

+ 0{\ 2 ) (h > t 2 ) . 



Although, in this case, the tree-level amplitude does not depend on the energy, we have 
nevertheless written it as A^\W) in eqs. (46) and (47), in order to be able to generalise 
the discussion to xPT in the following sections. 

In order to determine the K — > nn decay amplitudes we need to divide the correlation 
functions by the appropriate factors associated with the source (kaon) and the sink (two 
pions). For single particle states this is straightforward, and there are no finite volume 
power corrections, 

p -rn K \t K \ 

(0 \K 5 (0)KUt K )\ 0) = I/ 3 Z 2 K ^— , (48) 



2m 



K 



where the subscript indicates that the Fourier transform of the kaon field has been 
performed at zero momentum. 

The determination of the factors (0|7T7r|n) associated with the two-pion intermediate states 
\n) is more subtle. These matrix elements can be obtained by evaluating four-pion cor- 
relation functions. However, one has to be careful to project out the matrix element 
corresponding to the required state |n), which in the present case is a two-pion, cubically 
invariant state which contains an s-wave component. Although one can, in principle, try 
to use the behaviour of the correlation function with t 2 to project out the state with the 
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2m 



correct energy, this may be difficult in practice if this state has an energy close to one with 
no s-wave component. For this reason, it is useful to project out the contribution of the 

— # 

required state directly, by replacing exponential factors of the form exp(ik ■ x) in Fourier 
transforms by their s-wave projections, sin(Ax)/Ar (k = \k\ and r = \x\). For our purposes 
it will be sufficient instead, to average the correlation function C^-ft in eqs. (46) and (47) 
over all q\ 's which can be transformed into each other by elements of the cubic group (in 
the formulae below we actually choose to average over the v momenta with the same mod- 
ulus |<fi|). Such an average has a projection on the s-wave, but no projection onto 1 — 1,2 
and 3 components of the states. Since at one-loop order, both in the A<£> 4 theory and in 
XPT , I — 4 components (and higher partial waves) cannot be generated, such an average 
over cubically equivalent states is sufficient. We now carry out this procedure explicitly. 

We start by evaluating the averaged correlation function (see eq. (27)) 

Y^C^-^kMM) = EE (0kfc(*i)T-$i(*2)|n) ( n |CV(0)|m) (m \ K\t K ) | ) 

qi qi n , m 

^ e -™^l e -^p0|7r a (ti-t 2 )7r^(0)|W>(W|CV(0)|/O , ( 4 9) 

91 

where \W) represents the cubically invariant two-pion state with energy W which contains 
an s-wave component. As explained above, the sum over q ± , runs over all the three- 
momenta with the same modulus, \q\\. The ellipses in eq. (49) represent terms with a 
different behaviour in ti, and as always, we assume that we can use the time dependence 
to isolate the first term on the left hand side. 

In order to obtain the K — > ttti matrix element from the correlation function in eq. (49), 
we need to divide by the matrix element 

E(0 1^1-^)^(0) | WO . (50) 

91 

We can determine this matrix element from the four-pion correlation function: 
EE (°l T-ft(*2) ^(-t 2 ) nti-h) | ) 

91 V 

= e- 2Wt2 EE (° I - fe) 7r_ a (0) \W)(W \nl p (0) ^{-{t, - t 2 )) | 0) + • • • , (51) 

91 V 

where the sum over p is over the same range as the one over q\. Evaluating the four-pion 
correlation function in perturbation theory we find: 

E(0ka(*i)T-a(*2)^(-t 2 )7rjH 1 )|0> 

p,qi 



-2Wt 2 -2E( tl -t 2 ) g Xv , , 

V L y l ~ AEHJ) \ tl ^ t2 > 



and (52) 

(2E) 2 ^ V x 4£3 L 3, 

Combining the results in eq. (52) with those for the K — > tttt correlation functions in 
eqs. (46) and (47) we find that the one-loop contribution to the amplitude in both cases is 

(53) 



(W\Ow(0)\K) =Am(w) <_ x 



RefEim) + M2^EL + l?l^ 
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The finite volume corrections in eq. (53) are precisely those given by the Lellouch-Liischer 
factor (see subsec. 4.1.1). We stress that in the perturbative example studied in this 
section we use the conventional covariant normalization of states. The remaining diagrams 
in fig. 1 also contribute to the infinite volume result, but they do not have corrections 
which decrease as inverse powers of the volume. 

The above argument is a simple example of the discussion presented in ref. [5] , in which it 
was shown that from the correlation functions eq. (27), after dividing by the matrix ele- 
ments corresponding to the source and the sink, one obtains the modulus of the amplitude, 
up to finite-volume corrections given by the Lellouch-Liischer factor. The case £ 2 <C £i was 
analysed in ref. [4] where it was shown that the correlation function eq. (27) gives the real 
part of the decay amplitude up to power corrections in the volume which are not those 
described by the Lellouch-Liischer factor. This is before dividing by the two-pion matrix 
element eq. (50). Division by this matrix element removes a factor of cos(5) (where 5 is the 
s-wave phase-shift), turning the real part of the amplitude into its modulus, and modifies 
the finite- volume corrections into the Lellouch-Liischer factor [5] . The above exercise is an 
explicit one-loop demonstration of this effect, although at this order of perturbation theory 
the real part and modulus of the amplitude are equivalent (cos(5) = 1). 

The remaining diagrams in fig. 1 do not have singularities in the range of integration and 
hence only have exponential corrections in the volume (i.e. terms which decrease like, for 
example, exp(— mL)). To illustrate this feature consider the complete contribution from 
Ti in eq. (31), which is proportional to 



4.1.1 Evaluation of the Decay Amplitude 

The one-loop toy- model study in sec. 4.1 is an explicit illustration of the general procedure 
for evaluating infinite-volume amplitudes from correlation functions computed in a finite 
volume. Before proceeding to develop the discussion of sec. 4.1 in the context of chiral 
perturbation theory, we briefly summarise this general procedure for the determination of 
the decay amplitudes. For illustration we consider the case where the kaon is at rest and 
the total momentum of the two-pions is also zero (see however sec. 4.1.2). The procedure 
is as follows: 

1. We evaluate the correlation function C^-,^^, £i, £2) (defined in eq. (27)) averaged 
over all q\ with the same modulus. We envisage evaluating the matrix element 
with the two-pion energy equal to W = 2E qi + finite-volume corrections. Let 
Cq 1 -q 1 ,w(tK,ti,t2) be the component of C^-gi (tx, £1, £2) whose behaviour with £ 2 is 
given by exp(— Wt 2 ). 

2. We evaluate the four-pion correlation function (see eq. (51)) 




/ 



d 3 q 1 1 
(27r) 3 (2w) 2 2(w + E) 



+ 0(e-' mL ) . 



(54) 



^(£1, £2) = ( I TTftfa) 7T_ ?1 (£ 2 ) TTl - (-£2) 7T- 2 (-£l) | ) , 



(55) 
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where again we average over all q\ and q 2 with modulus equal to |<fi|. As above, let 
C^ 2 W (ti,t 2 ) be the component of C^ 2 (t 1: t 2 ) whose behaviour with t 2 is given by 
exp('-2W* 2 ). 

3. Finally we calculate the kaon propagator, C^{t K ) 

C$(t K ) = (0\K s (-t K )Kl(t K ) |0> (56) 

for sufficiently large \t K \ so that the correlation function is dominated by the kaon- 
state. 

4. The finite-volume matrix elements are given by: 

| v (nn\O w (0)\K)v \ = C *-*MtK,h,t 2 ) = (5?) 

^c^ w {t u t 2 )q{t K ) 

where the subscripts V on the matrix element denotes that it is the finite-volume 
matrix element (with the finite- volume normalization). Note that the K — > im 
matrix elements discussed in sec. 4.1 in general and eq. (53) in particular, which 
were presented using the covariant normalization, should be divided by \ / 8i71kE 2 V 3 
to translate them into the finite- volume normalization in eq. (57). 

5. To obtain the infinite-volume matrix elements {7nr\Ow(ty\K) (with the conventional 
covariant normalization) from the finite-volume ones, we multiply the latter by the 
LL factor [3, 4] 3 

\(nn\O w (0)\K)\ 2 = W 2 {q<P'{q) + kS'(k)} (^f \ v (7nr\O w (0)\K) v | 2 , (58) 



where k = yW 2 /4 — m 2 , q = kL/(2ir), 5(h) is the s-wave (and in our case 1 = 2) 
ix-ix phase-shift and the geometrical function (p(q) is given by 

2ix 2 a 

ta " W ^ - S az a/ ( rw) ■ (59) 

This procedure holds for the matrix elements discussed in this paper in full QCD and also 
in the quenched approximation (this is not the case in the quenched approximation for 
AI = 1/2 transitions [8]). The example in sec. 4.1 is a model one-loop demonstration of 
the validity of the procedure. 

At present the LL factor has not been generalized to the situation in which the two-pion 
final state has a non-zero total three-momentum, and so we are unable to carry out step 
5. We now discuss this point in a little more detail. 



3 In K + — > 7T+7T decays the final state particles are distinguishable and so the factor in eq. (58) differs by 
a factor of 2 compared with the corresponding one presented for indistinguishable final state particles [3,4]. 
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4.1.2 Finite- Volume Corrections in a Moving Frame 

In order to implement our strategy, which requires the computation of matrix elements 
at unphysical kinematics in general and SPQR kinematics in particular, simulations have 
to be performed with the two pions having non-zero total three-momentum. Thus, if the 
calculations are to be precise up to exponential corrections in the volume, we need the 
generalization of the Luscher quantization condition and the LL factor away from the 
centre-of-mass frame. The quantization condition in a moving frame has been presented 
in ref. [36], whereas the generalization of the LL-factor is being studied and we will report 
the results in a future publication. We anticipate that the generalization of the LL-factor 
will be obtained from the quantization condition in an analogous way to the derivation in 
the centre-of-mass frame given in ref. [4]. 

At this point we are able to check whether the energy shift obtained in perturbation theory 
following the procedure of sec. 4.1 generalized to a moving frame, agrees with the expansion 
of the quantization condition presented in ref. [36]. This is indeed the case. Repeating steps 
(34) - (39) for two particles with energies E\ and E 2 we find that the energy shift is equal 
to \v d I (AE1E2L 3 ), where v d is the number of terms in the momentum sum with energy 
equal to E\ + E 2 . This agrees with the result obtained by expanding the quantization 
condition in ref. [36] and using the one-loop expression for the phase-shift 4 . 

We postpone the general discussion of finite-volume corrections at non-zero total momen- 
tum until the corresponding LL-factors have been derived and now proceed to study one- 
loop chiral perturbation theory. 



4.2 K — > 7T7T Decays in on a Finite Volume 

In this section we generalise the previous discussion to %PT . In this case the weak am- 
plitude in general depends on the energy, and hence other finite-volume correction terms 
appear in the correlation function. To illustrate this, consider the term in eq. (31) which 
contains T 2 , but now generalized to allow for the energy dependence in the lowest order 
weak and strong amplitudes (the power corrections arise from singular terms as explained 
in the previous section): 

{1 2(E-w)t 2 _ I 1 

WE) - j- 3 £ 2{E _ w) A(2w)Af(w) I (60) 

where A(2E) is the tree level amplitude and M(w) is the strong interaction amplitude (in 
the model of sec. 4.1 J\f(w) = \/Aw 2 ). Separating the contribution with w = E from the 
rest, we find 

1 e 2(E-w)t 2 _ 1 

k-.w^E v ' 
We thank Kari Rummukainen for helping us to clarify this point. 



C = g(°\t K ,t 1 ,t 2 ,2E){A(2E) 



uN{E) 
IT" 
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g<-°Ht K ,t u t 2 ,w) 



W E ) - 73 £ 



g[w) 



L 3 ^ 2(w 2 - E 2 ) 



(61) 



where 



g (w) = (w + E)A(2w)Af(w) , (62) 

and the ellipses represent terms which have a different t 2 behaviour (as above, we assume 
that we can isolate the terms with the required behaviour in t 2 ). 

We now use the summation formula (for smooth functions F), 



1 v / F(k 2 

L 3 ^ K 2 

k 



k' 2 



7 



d 3 k F{k 2 



(2vr) 3 K 2 - k 2 



(63) 



where the prime on the summation indicates the omission of the v terms with k 2 = K 2 and 
F'(K 2 ) denotes the derivative of F with respect to k 2 , evaluated at k 2 = K 2 . Eq. (42) is an 
example of this summation formula with F(k 2 ) = —l/(4w) and w 2 = k 2 + m 2 . Applying 
eq. (63) to eq. (61) we find 



C = g^ \t K ,t 1 ,t 2 ,W)S [ A(2E)(l-h(E)) + 1 ^^^- 



w=E 



-Q«>\t K MM, W)A(2E) {EN(E)^- L - ^ ^ + EW(u,) } 



dw 



, (64) 



w=E 



where —h(E)A(2E) is the one-loop correction to the infinite-volume amplitude with 

d 3 k g{w) 



HE) = \vj 



(2tt) 3 w 2 - E 2 



(65) 



The one-loop results for the AI = 3/2 matrix elements of 04,7,8 are presented in detail in 
Appendices B and C and on the web-site [23]. We now look at the finite- volume corrections 
on the right-hand side of eq. (64). The term proportional to dA/dw has the role of shifting 
the argument of the amplitude from 2E to W = 2E + uAf(E)/L 3 . In an earlier study this 
term had been misinterpreted as a genuine finite-volume correction, whereas it is present 
to shift the argument of A. For this reason the authors of ref. [27] concluded that there are 
different finite-volume corrections for the two cases when rnx = rn-K and uik = 2m n , with 
the two pions at rest. These corrections must instead be the same since they only depend 
on the final state. Our explicit calculations, presented in sec. 4.3 below, demonstrate that 
this is indeed the case. 

The second line on the right-hand side of eq. (64) is a universal correction which depends 
only on strong interactions and gives the Lellouch-Luscher factor. This concludes our 
demonstration of finite-volume effects in one-loop %PT . 
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4.3 K + -> tt+tt Decays in xPT at Finite Volume 

We now briefly summarize the results obtained in full QCD in finite volume for the physical 
amplitudes, i.e. for a kaon at rest decaying into two pions with the insertion of a weak 
operator at zero four-momentum transfer. Of course, our approach requires us to perform 
lattice simulations (quenched and unquenched), with the insertion of non-zero momentum 
at the weak operator and the corresponding results will be presented below in sec. 5.2. 

For the calculation of the correlation functions needed to determine the K + — > -k + tt° decay 
amplitude at one-loop in x?T we need the following propagators: 

g W^\t X ty | 



/ 



Jd'xe^-^i^Q^^ty)) = 6 ^ - (66) 

i e -Wl\tx-ty\ 



where w| = m 2 + \q\ 2 , e represents the sign function, and the index j runs over the spatial 
components 1,2 and 3. 

For the evaluation of the tadpole diagrams we need the relations: 

d 4 k 1 1 1 

— r > and 



(2tt) 4 k 2 + m 2 L 3 - 2w ? - 



(i 4 A; fc 2 1 ^ / r dk A m 2 \ 

(2tt) 4 k 2 + m 2 = 2^ _ 2^J' ( ^ 

which are valid up to terms which are exponentially small in the volume 5 . 

When evaluating the diagrams (6), (c) and (d) we follow similar steps to those in eq. (31). 
In xPT a number of integrals appear with different numerators. The complete set of these 
integrals and the corresponding results, obtained upon integration over A; 4 and / 4 , are given 
in Appendix A. 

The result for the matrix element of O4, extracted from the correlation function with the 
kaon at rest and with W = m K , where W is the finite- volume two-pion energy, is given by 
the expression: 



(7T + 7r O |04|^ + > 



6V2a^ 27 ^ 



f K f 2 

2 (l a + I b E + I c+d ) 



4 

m K 



(47T) 2 / 2 



(4tt/)= 

+ 0(m A KlT ) counterterms , (68) 



where y/s is the invariant mass of the two-pion final state. In this case \fs = uik, however 
since the finite-volume corrections depend only on the energy of the two-pion state, we 

5 The divergent terms proportional to J dfc 4 /(27r) cancel when all contributions to the amplitude are 
summed. 
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write the answer in terms of y/s. The functions I z f t a,b,c+d are defined in eqs. (80) and (81) 
below and 1^ = ReI b . a^ 27 ' 1 ) is the leading order low-energy constant, and new constants 
appear at order m\ K . As explained in the introduction, in our approach it is envisaged that 
all the low energy constants up this order will be determined by studying the behaviour of 
the matrix element as a function of masses and momenta. F(\fsL) is the centre-of-mass 
finite volume correction 

The shift in the energy of the two-pion state due to the finite volume is given by 

AW = W-2E= v 8 1 . (70) 

From eq. (70) with the pions at rest, using AW = — Ajra 1 ^ 2 / (m n L 3 ) [34], we can obtain 
the infinite- volume I = 2, n — n scattering length 

1=2 m n /"-i \ 

°° =~^P- (71) 

We remain now with the limit rriK — > 2m n , corresponding to the two pions at rest, and 
compare our results with those of ref. [27]. In this case 2z n — > 2z = —17.82726584 and 
v — 1 so that the finite volume correction becomes 

This result is in agreement with the general LL formula, whereas it disagrees with the 
result in ref. [27] (see eq.(3.3) in [27]), where the corresponding correction was found to be: 

17.827 5tt 2 

4m n L + 2m3L 3 ' ^ ' 

We have explicitly verified that the discrepancy is due to the fact that in ref. [27] it was 
not recognized that part of the 1/L 3 finite- volume correction is proportional to dA/dw 
and is absorbed by the shift of the argument of the amplitude, A(2E) — > ^4(W) (see the 
discussion in section 4.2 above). This is also the origin for the erroneous conclusion of these 
authors that the finite- volume effects are different for the two-pions at rest with mx = 2m n 
and rriK = m n [26, 27]. 

For the electroweak penguin operators 7 and O s we obtain 



(7r+7r°|07,8l^~ 



1 + ^fy 2 d,f + FWiL)) + j0j- 2 (J! + J! + J? +d ) 
+ 0(m 2 K 7r ) counterterms , (74) 



23 



where in this case y/s = tuk, 7^ 8 ' 8 ^ are the leading-order low-energy constants (for Oj and 
O s ) and the functions Jf are the real part of the functions Jj defined in eq. (87). As 
expected, the finite volume corrections encoded in F(y^sL) are the same as for O4. 

For K — > mi matrix elements for which the two pions have zero total momentum, the 
finite-volume corrections take the universal form 

(n + n°\O w \K + ) = LO + jJL^F(y/5L) + • • , (75) 

where LO is the lowest order contribution and the ellipses represent NLO infinite-volume 
terms. For AI = 3/2 transitions this is also true for the quenched approximation, since 
the relevant one-loop contributions are as above. For AI = 1/2 transitions this is not the 
case [7, 8]. In the following section and in appendices B and C we present the results also 
for unphysical kinematics at one-loop order in x?T without explicitly exhibiting the finite- 
volume effects. As discussed in sec. 4.1.2 the general theory of finite-volume corrections 
to matrix elements with the two-pions at non-zero total momentum has not yet been 
developed. We are currently investigating this question, and will present the perturbative 
finite-volume corrections within the general framework. 



5 One-loop corrections in 

In the preceding section, we have demonstrated that, at one-loop order in chiral pertur- 
bation theory, the finite-volume corrections to the amplitudes extracted from correlators 
in which the pions are annihilated at the same time or at different times are exactly the 
same and correspond to the general LL formula. As explicit examples, we have given the 
one-loop amplitudes for the AI = 3/2 operators, O4, 7 and C 8 , in the physical case (i.e. 
with no momentum insertion at the weak operator), for a kaon at rest decaying into two 
pions in full QCD. In order to implement the strategy described in sec. 2 however, we 
also need the amplitudes in other kinematical situations for which there is an insertion of 
momentum at the weak operator. In this section we discuss the evaluation of the one-loop 
corrections in the chiral expansion, both in full QCD and in the quenched approximation 
for arbitrary meson masses and momenta. The general results are very long and compli- 
cated: for full QCD they are presented in appendix B.2 (and also on the web site [23]) 
and for the quenched approximation on the web site [23]. In presenting our results we do 
not exhibit explicitly all the NLO counterterms, i.e. the terms proportional to the /3j's in 
sec. 2.1 and the <Vs in sec. 2.2. Of course these terms must be included in any comparison 
of lattice data with %PT . 

We now consider separately the following special cases, which are useful for current nu- 
merical simulations: 

1. physical kinematics with the kaon at rest and with the mesons having arbitrary masses 
rriK and m^. By "physical" here we mean that there is no insertion of momentum at 
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the weak vertex, so that the two pions each move with momentum k = ym|/4 — vn\. 
The results for this case are presented fully in sec. 5.1 below; 

2. SPQR kinematics with the kaon and one of the final-state pions at rest and the other 
one moving with arbitrary energy E n . The final state is symmetrized to ensure that it 
has 1 = 2. The results are presented in appendix B.l for full QCD and appendix C.l 
for the quenched theory; 

3. choices of the quark masses such that tuk = m n and tuk = 2rn n , in both cases with 
all the particles at rest. These are presented in sections 5.2.1 - 5.2.4. 

Below we follow the notation introduced in the preceding sections and assume SU (2) isospin 
symmetry (m s ^ m u = ma), together with the following relation for the rj mass: 

m l = \™\ - \ m l • (76) 

Ultra-violet divergences are regulated using dimensional regularization and, following the 
conventions of ref . [20] , renormalized by performing the subtraction of the term 

l - - lE + log(47r) + 1 . (77) 



5.1 Matrix Elements for the Decay K + — > 7r + 7r° with Physical 
Kinematics 

In this section we describe the calculation of the matrix elements for the decay K + — > 7r + 7r°, 
in both full QCD (sec. 5.1.1) and the quenched theory (sec. 5.1.2). We underline again 
that by physical we mean that there is no insertion of momentum at the weak operator. 
The masses of the mesons are arbitrary. 



5.1.1 Calculations in Full QCD 

We start with the discussion in Minkowski space. At lowest order the matrix element of 
04, obtained using the chiral operator defined in eq. (12), is given by 

= - 7^ « - ml) . (78) 

At one-loop order we have to compute the diagrams of fig. 1. In the tadpole diagram (a), 
all pseudoscalar mesons propagate in the loop. In the remaining diagrams, because of the 
symmetries of the Lagrangian, it is the combinations indicated in tab. 1 which contribute. 

Our results for this case agree with those of ref. [37] and with those of Golterman and 
Leung [26], except for the finite terms which depend on the regularization/renormalization 
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Diagram 


Particles 


(b) 


7T + 7T 


(c) 




(d) 


KV, K+1T-, K°T] 



Table 1: Combinations of mesons looping in the different diagrams. 



(and the sign of the imaginary part of diagram (b) in fig. 1). The matrix element is given 
by 



(7T + 7r |C 4 |^ 



tefl 



[m\ - ml) 



1 + 



m 



(4tt/) 



K j 



+ 



m 



K 



(4tt) 2 / : 



(I a + I b + I c+d ) 



(79) 

plus the 0{p A ) counterterms (i.e. the terms proportional to the /3j's in eq. (6)). The 
integrals in eq. (79) are given by 
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V18y 2 9y ' 18 3 ; 
where y = m\fm 2 K , \x is the renormalization scale, 

A( y) = log (£$=lf ) and ' 

'7 - 4y + VI -44y+ 16y 2S 



(81) 



1 _ 44y + I6y 2 log 



= -2J44y - 1 - 16y 2 arctan 



7 - 4y - VI -44y+ 16y 2 / 
/ V44y-l-16y 2 \ 

V 7-42/ 



(82) 



In deriving the above expressions we have used eq. (76) and the fact that, with physical 
meson masses, the argument of the square root in eq. (82), 1 — AAy + 16y 2 , is negative. 
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From eq. (79), we can extract the phase of the amplitude, 8, defined by 

A=\A\e iS . 

We find 



5(k = \Jm 2 K /4: — to 2 ) 



m 



K 



;i-2yWl-4y 



(83) 



(84) 



16tt/ 2 

Rewriting the right-hand side of eq. (84) in terms of k we recover the expression for the 
1 = 2 strong interaction n-n scattering phase, 



S(k) 



(2k 2 + ml) k 



8nf 2 E ' 

as required by Watson's theorem. The 1 = 2 scattering length is given by 



(85) 
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= lim 



6(h) 



k 87r/ 2 ' 

For the electroweak penguin operators O-j and 0$ the corresponding formulae are: 

2 V / 2 7 (8,^ 



(86) 
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(87) 



where I z f is given in eq. (80) plus the 0(p 2 ) counterterms (i.e. the terms proportional to 
the S^s in eq. (9)). The integrals appearing on the right-hand side of eq. (87) are 
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(88) 



^4(y) and -B(y) have been defined in eq. (82). The strong interaction phase is of course, 
the same as for O4, since it only depends on the isospin of the two-pion state. For this 
case our results agree with those of ref. [38]. In ref. [24] the one-loop contributions to the 
matrix elements are presented numerically and we agree with these results. 
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Calculations in Euclidean Space The calculation of the Euclidean matrix elements 
in terms of the low-energy constants at NLO in %PT proceeds along similar lines. In the 
Euclidean theory the imaginary part is zero and the arguments of the logarithms correspond 
to the absolute values of the arguments in the Minkowski case. For a generic (Euclidean) 
kaon momentum px (which is also therefore the momentum of the two-pion state) the 
logarithms take the form 

/ 



log 



V 





A+4f) 


\ 




/l + 4f) 





Setting p 2 K — > —m\ we recover the Minkowski amplitude without the phase 
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The functions I x are those defined in eqs. (80) and (81) and Jjf = Re I b . 



I a + 7f + I c+d ) 
(90) 



5.1.2 Calculations in Quenched QCD 



In this section we present the results for the matrix elements in quenched QCD. As above 
we do not exhibit the NLO counterterms given in eqs. (6) and (9) explicitly but of course 
they must be included in any comparison with lattice data. 

The matrix element of the operator O4 takes the form: 
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with 



vy; v y y & v 3 - 2 y- VI -8y + 4y 2 ; 



For the EWP operators the corresponding results are 
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where 1^ is defined in eq. (92) and the remaining integrals are given by 



(94) 



(95) 



7 q 



J. 



c+d 



+ 



i-y 



m 



+- 



a 



logl -f I + -log(y) + 



9m 2 



A 



-2 + -Ugf^ 



i-y 



y 



1 - ^y ) log! 



y 



m; 



4 , y 



3 3 



m 



y 



+ (1 - 2y) 



- m Jl -Ay -I 



1 1 



= - + £ log - 1 - -A(y) - - + - log(v) 



771 



A 



^(y) 



2y 



2y, 



■l + 6y- 16y 2 + 8|/ 3 



12y 2 (l-y) 12y 2 (l-y)(l-82/ + 42/ 2 ) 

+ 



C(y) 



4 1 , 3-6y-2y 2 1 (2-y 

+ ~ ~ — + "TT^ nVZ -log 



+ a 



9 3y 3Qy 2 (l-y) 
1 -Qy + 4y 2 ^ 



12y 2 (l-y) 



^y) + 



v y / 
i i 

%" 9 



17 



+ U + V~y~ % + ^ 



c(y) 



(l_ y )(l_ 8y + 4y 2) 



+ 



+ 



29 



9 - 3Qy + 42y 2 - 28y 3 + 8y A (2 - y 
36y 2 {l-y) ° g \ y 



(l-y)(l-%) 
3y 2 



log(y) 



(96) 



5.2 K + -> tt+tt Matrix Elements with the SPQR Kinematics 

As explained in the Introduction and section 2, our aim is to extract the physical K — > 7r7r 
matrix elements from simulations at unphysical kinematics using NLO . This requires 
the calculation of the one-loop chiral corrections to the K + — > tt + tt° amplitudes with the 
insertion of weak operators carrying four-momentum. In this section we consider the matrix 
elements for the SPQR kinematics, in which the kaon and one of the final-state pions are at 
rest and the second pion has an energy E n . We have seen in section 2 that the computation 
of these matrix elements are sufficient to determine all the required low-energy constants 
at NLO in the chiral expansion. 

In order to ensure that the final state has I — 2, it is necessary to symmetrise over the 
momenta of the two pions. Note also that the only diagram which contributes to the 
rescattering phase is diagram (6). 

The results for full QCD are presented in appendix B.l and for the quenched theory in 
appendix C.l. Here we present the results for the two choices of quark masses corresponding 
to m K = 2m n or m K = m w in each case with both pions at rest. 



5.2.1 Matrix Elements with m K = 2771^ and both Pions at Rest in Full QCD 



Consider first the operator C 4 . In the limit where both pions are at rest (cu — 1) and 
rriK = 2m n (z = 1/2) from eqs.(108) - (111) we can readily recover the corresponding 
result (y = 1/4) in eqs. (79) - (81): 
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The corresponding result for the EWP operators is: 
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5.2.2 Matrix Elements in the SU(3) limit, m K = m„, and with both Pions at 
Rest in Full QCD 

In the SU(3) limit m s = m d = m u , for which m K = m n , with the two pions both at rest 
eqs.(108) - (111) reduce to 
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in agreement with the result given in eq. (82) of ref. [26], up to scheme-dependent finite 
terms. 



For 07 i8 the corresponding result is: 
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5.2.3 Matrix Elements with m K = 2m 7r and both Pions at Rest in Quenched 
QCD 



The results for the matrix elements with the SPQR kinematics in the quenched theory are 
presented in detail in appendix C.l. In this and the following subsection we present the 
results for the special cases m K = 2m w and m K = m n respectively, in each case with both 
the pions at rest. 

For the matrix element of in the quenched theory with rnx = 2m n We find 
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We are grateful to M. Golterman for revisiting this calculation and confirming his agreement 
with the expression in eq. (101) [39] (the result in eq. (101) is in disagreement with that 
presented in eqs. (3.1) and (3.2) of ref. [27]). The corresponding result for the EWP 
operators is: 
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5.2.4 Matrix Elements in the SU(3) limit, m K = m^, and with both Pions at 
Rest in Quenched QCD 



Finally we present the results with rriK = m w m the quenched theory with the pions at 
rest. In this case for the matrix element of O4 we do agree with the result in ref. [26] 
finding 
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For the matrix elements of O-j^ the corresponding result is: 
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6 Conclusions 

The precise evaluation of K — > im matrix elements in lattice simulations is a major long- 
term challenge. In this paper we have explained our strategy for the reduction of an 
important source of systematic error in current studies, viz. the use of xPT at leading order 
in which final state interactions are neglected. We propose to use lattice computations of 
K — > im matrix elements with meson masses and momenta which are unphysical (but 
accessible to lattice simulations) to determine all the low energy constants necessary to 
obtain the physical amplitudes at NLO in x?T . We have presented all the necessary 
ingredients to implement this strategy for AI = 3/2 transitions. In particular we have 
evaluated all the chiral logarithms at one-loop order for arbitrary masses and momenta. 

We have shown in sec. 2 that it is possible to determine all the necessary low-energy 
constants from matrix elements with the SPQR kinematics, i.e. with the kaon and one of 
the pions at rest and with the second pion having an arbitrary momentum. We present 
the results for this choice of kinematics separately, but stress that of course this is not the 
only choice for the determination of the low-energy constants (for this reason we present 
the very lengthy results for general kinematics). 

We incorporate into our work the recent progress in understanding finite-volume effects 
in K — > Tire decays. For AI = 3/2 transitions we find that these effects are given by the 
Lellouch-Luscher factor (for zero total three- momentum), even in the quenched approx- 
imation. This is not the case for AI = 1/2 transitions, where the absence of unitarity 
in the quenched theory leads to major subtleties and difficulties, such as the dependence 
of the final state interactions on the choice of AI = 1/2 weak operator and non-standard 
behaviour of the correlation functions with time and volume. This is the subject of a paper 
under preparation. We are also investigating the generalization of the LL-factor to the case 
in which the total three-momentum of the two-pions is not equal to zero. 

In ref. [17] a particular implementation of our proposed strategy is discussed. The conclu- 
sion of ref. [17] is that it is possible to determine all the low-energy constants necessary to 
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obtain the physical matrix elements for (27,1) and (8,1) operators from the knowledge of: 



i) K — > 7T7T matrix elements with both the pions at rest for nix = 2m n and nix = m n ; 

ii) K — > n matrix elements at non-zero momentum; 
iii) the matrix elements for K°-K° mixing. 

We agree with this conclusion for the AI = 3/2, (27,1) operators studied above. Notice 
however that our complete basis for the %PT weak operators at order p 4 differs in the 
replacement of the operator O20 with our operator O22 ■ The two choices turn out to be 
equivalent also in the most general case with unphysical kinematics (i.e. with an injection 
of momentum by the weak operator). 

For AI = 3/2 decays our strategy can also be implemented for partially quenched lattice 
QCD, and we will present the corresponding one-loop results in a future paper. (For the 
case rriK = m^, with all the mesons at rest and with degenerate sea-quarks, the matrix 
element of C 4 has been evaluated in partially quenched QCD in ref. [28].) 
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A Integrals which appear in the evaluation of the cor- 
relation functions 

In this appendix we list the energy integrals which appear in the evaluation of the contri- 
butions of diagrams (b), (c) and (d) of Fig.l and present the corresponding results. We 
write each integral in the generic form: 




(105) 



6 The numbering of operators follow ref. [22] in both cases. 
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where Af is a polynomial in fc 4 and U, w\ = m\ + k 2 and w\ = m 2 , + P. 
Using the identity 

r dx x n e~ txt /-^V /* ^ e_lx * 

7 2^x 2 + w 2 = \di) J 2^x 2 + w 2 ( ' 

we obtain 



V / 


_ 1 c -(t«l+t«2)|t| 

(2 Wl )(2w 2 ) 




„• ^ „— (u>i+W2)\t\ 

= —i£(t) e v ;M 

4w 2 


J(/ 4 ) 


= —iE\t) e v ;M 

v 7 Awi 


/(M 4 ) 


1 -(tui+iu 2 )|t| 

4 




_ 5 W w l c -(»l+» 2 )ltl 

2w; 2 4w 2 




2wi 4u>i 




= , 5 (t)^ e -(-i+- 2 )l*l 
4 


/(M 4 2 ) 


4 


mil) 


= m J 2n 2 



:tt 2 ° l j+ 4 

where e represents the sign function. 



B Results in Full QCD 

In this appendix we present the results for matrix elements in full QCD. The results for 
general kinematics, i.e. for arbitrary quark masses and momenta, with the two pions in 
an 1=2 state, are given in subsection B.2. We start however, with subsection B.l which 
contains the results for SPQR kinematics. 

The integral I z f, which corresponds to the one- loop contributions from the renormalization 
of the mesonic wave-functions and from the replacement of the factor 1/ f 3 which appears 
at lowest-order in the chiral expansion, by l/(fxf%) is given by (y = m 2 /m 2 K ): 

/ m 2 \ 7 1 fm 2 \ 

I zf = -3(1+ y) log + - log(y) + - (y - 4) log [-jj . (107) 

The results are given in terms of the variables, y = m 2 /m 2 K , z = m^jm^ and u = E n /m n , 
where E n is the energy of the pion whose momentum is (in general) not equal to zero. 
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B.l K+ -> tt+tt Decays with SPQR Kinematics in Full QCD 

B.l.l Results for the matrix elements of 4 in full QCD with SPQR kinemat- 
ics 

The matrix elements of O4 with SPQR kinematics is given by: 



{7r+n o mK+) = Z 6v^ (27>1) f g (l+a; + 2z) / 
V ' ' ' fxf 2 \ 2cu \ 



m K 



16tt 2 / 2 

where I z f is given in eq. (107) and I a ,b,c+d are as follows: 



16tt 2 / 2 

(I a + I b + I c+d ) ) , (108) 



G 1 



1 4 (l2 z 2 - uj + 8 z (1 + uj) \og(z 2 



- (5Qz 2 -Alo + Uz (1 + uj) + 2Qz 3 (1 + uj) +4z 4 (21+wj) log 

+ z (-4 + ^ 2 ) (3 + 2z + 3uj) log ^^§j |, (109) 



3(l + 2z + w) + 3ji-p^(l + 2z + w)log _ n — 

\ 1 + A/ l +a j 



+ 



,2 ' 



-3 - 2cj + uj 2 - 2z(3 - 3cj + uj 2 )] log f ^ 1 | , (110) 



and 



Ic +d = Il +d + ll + Aog (^j + Il +d \0g{z 2 ) + It +d log + Il +d . (Ill) 

The separate components in eq. (Ill) are as follows: 

1 2 z 5 to 2 (-16 - 15 uj + uj 2 ) + 3uj 2 (3 - 26 uj + 3 cj 2 ) 



/ z 2 



c+d ' I -2(_0 V J- ,., J- *2. ■ ^ 2 



9uj 2 (-2z + UJ + Z 2 UJ 
+2zlu (-18 + 140cu-33w 2 + cj 3 ) + z 4 uj (43 + 129 uj - uj 2 + 9 uj 3 ) (112) 

+z 2 (36 - 269 uj + 174 w 2 - 79 uj 3 + 18 u; 4 ) - 2^ 3 (ll + 78 uj - 107 uj 2 + 48 uj 3 + 30 cj 4 ) J , 

P c+d = ^—^ I -2uj 2 -2z 4 uj 2 + 4zuj (1 + uj) + z 2 (-6 + 28 uj- 6 uj 2 ) 

(-l+oj 2 )], (113) 



-z 3 
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\54 (-! + *) uj 2 {-2 z + uj + z 2 ujf , 



,uj 2 



| - 18uj 5 + 18zuj 4 (7 + 2uj) +Az w uj 3 (-7-3cj + 4( 

-6 z 2 uj 3 (63 - 11 uj + 21 uj 2 ) - 2 z 9 uj 2 (-135 - 50 uj + 12 w 2 + 8 uj 3 ) 
+5z 3 uj 2 (126 - 179uj+ 129uj 2 + 8uj 3 ) + 3 z s uj (-90 - 156u - 265uj 2 + 71 uj 3 ) 

+z 4 uj (-576 + 1926 uj - 1373 cj 2 + 75 uj 3 - 202 uj 4 ) (114) 
+/ (36 + 582 + 1116 uj 2 - 357 uj 3 + 1047 u; 4 - 36 uj 5 ) 
+6 ^ 5 (36 - 179 uj + 222 uj 2 - 210 cj 3 + 197 uj 4 + 56 cj 5 ) 

-6z 6 (42 + 85cj- 195cj 2 + 309cj 3 + 29c/ + 77cj 5 ) J, 



I 4 



r+d """" U08 (-1 + *) uj 2 (-2 z + uj + z 2 ujf J 

x | - 36 z 10 uj 4 - 36 uj 4 (1 + uj) + z 9 uj 3 (l87 + 3uj - 4uj 2 

+12zuj 3 (l8 + 29 uj + 3 uj 2 ) + z s uj 2 (-378 + 11 uj + 141 uj 2 + 4cj 3 ) 
-r 2 w 2 (432 + 1052 uj + 303 cj 2 + 163 uj 3 ) 

-3z 7 uj (-114 + 6cj + 188cj 2 + 29cu 3 + 21cj 4 ) (115) 
-3 z 4 uj (230 - 30 uj + 363 uj 2 + 151 uj 3 + 28 uj 4 ) 
+z 3 uj (288 + 1314 uj + 566 w 2 + 771 uj 3 + 163 uj 4 ) 
+3 z 6 (-36 - 26 uj + 288 cj 2 + 14 uj 3 + 89 uj 4 + 21 w 5 ) 

+3 z 5 (36 - 178 + 192 uj 2 - 111 uj 3 + 19 + 28 uj 5 ) J 



i 



6w 2 (-2^ + cj + 2 2 w ) 3 / Vi + v/T 



a; 



2 



X 



|l8 2 5 c; 2 + 6cj 2 (-1 +4cj) -3^ 4 cj (5 + 2cj + lieu 2 ) -3^cj (-8 + 21cj + 5cj' 
+3^ 2 (-8 + 17 uj - 4uj 2 + 7 uj 3 ) +z 3 (-2 + 24^-11^ + 25^ 



2z 2 (1 + z) 



+ I o7 ) V(-! + *)* (-2 - * + z 2 ) x (116) 

\9 (-1 + z) uj J 

, /2 + 3^-2^ 2 -2a/(-1 + 2) 2 (-2-z + z 2 )\ 

log V + 

\2 + 3z-2z 2 + 2V(-l + z) 2 (-2-Z + 2 2 )/ 
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12 zoo- 12z 3 u-8 uj 2 + 4 z 4 uj 2 + z 2 (9-5 uj 2 ) x 



{' 



v 54cj 2 (-2^ + cj + 2 2 cj) 3 / 

6uj 3 + 3z 4 uj (-11 + 5cj 2 ) +zu 2 (-31 + 13cj 2 ) + z 2 (33cj-21cj 3 ) 

+z 5 (14 cj 2 - 8 uj 4 ) + z 3 (18 - 19 uj 2 + 13 cj 4 ) j 
'3 z + 2w-2 z 2 cj - ^12^^- 122 3 c;-8cj 2 + 4^ 4 c«j 2 + 2 2 (9-5cj 2 



log 



?,z + 2uj-2z 2 u] + J12zuj- 12 z 3 uj -8uj 2 + 4z a uj 2 + z 2 (9-5uj 2 ) 



B.1.2 Results for the matrix elements of (9 78 in full QCD with SPQR kine- 
matics 



The matrix elements of the EWP operators are given by: 



(«+«°\O rfi \K+) = |^ (8 ' 8) 



1 + 



m K 

16vr 2 f 



(I z f + J a + J b + J c+d ) 



(117) 



where I z f, the one-loop contribution from the renormalization of the mesons' wavefunctions 
and the replacement of 1/ f 3 by 1/ fxfn * s given in eq. (107) and the J a ,b,c+d are as follows: 

J a = \ |l8 log(, 2 ) -2 (11 + 10 z 2 ) log (^pj + (-4 + z 2 ) log (^pj } (118) 

(119) 



l+Ul 



and 



J, 



-(1+^(^ + 4^-^(5 + 3^) ( -2uj-Qz 2 uj+15z(1+uj) 

/ ^ o \ "T"( ^ 



+ 



c<j(— 2,2 + + z 2 u;) 
1 



6cj J \ //- 

2cu 3 -4 2 6 cj 2 (1 + cj) + 3^c; 2 (3 + 5w) + 



6(-l + ^(-2,2 + a; + z 2 lj) 2 / 
2z 5 uj (l5 + 10cj + 7cj 2 ) -z 2 cj (34 + 77cj + 13 cu 2 ) - z A (30 + 46 w + 93 uj 2 + 5 cu 3 ) 

+ z 3 (30 + 102 cj + 41 cu 2 + 43 uj 3 ) J log(^ 2 ) + 

N ]8cu 3 -4/cj 3 -42cj 2 (7 + cj) + /cj 2 (17 + 5cj) + 



,12 (-1 + z) uj (-2 z + w + z 2 uj) 

z 5 uj (-28 - 21 uj + uj 2 ) +2 z 2 uj (l6 + 6 + 5 uj 2 ) - z 3 (l8 + 16 uj + 25 cj 2 + cj 3 ) + 

(l8 + 32 + 5 uj 2 + 5 uj 3 ) | log ^^fj 



(120) 
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'Quj + Qz 2 lu-z (7 + 5u 2 )\ r .I fl-VT^ 



2u{-2z + u + z 2 uf ) {l + VT^J 2 



2£ \ [—— ^TT- -17, /2 + 3z-2z 2 ^/(-l + z) 2 (-2-z + z* 

2 ) V (- 1 + ^) (- 2 - * + * 2 ) lo S 



>3(-l + z) j ' ^ 2 + 3^-2^ + 2 v /(-l + ^) 2 (-2-2 + 2 2 

+z f 2 ^ + ^^~^ (3 + ^ 2 2) l x ^12^-12,3^-8^ + 4^^2 + ^2 (9-51,2) 
\ 6uj (— 2 z + + 2 2 (J) J 

'3z + 2u-2z 2 u - yjl2zu- 12z 3 u - 8u 2 + 4z 4 u 2 + z 2 (9-5u 2 ) 



x log 



3z + 2lu-2z 2 lu + ^12 zuo- 12 z 3 lu - 8lu 2 + 4 z 4 lu 2 + z 2 (9-5cj 2 



B.2 Results for the Matrix Elements in Full QCD for General 
Kinematics 

In this section we present the results for the matrix elements in full QCD for general 
kinematics, i.e. with arbitrary masses and momenta for the mesons. The results are rather 
long and complicated so we also present them on the web site [23] . The results are presented 
in terms of the following functions: 

1 

72 ' 



dip) 



P 

div(p) = -1-log^ 



B(p,mi,m 2 ) = y (1 — m 2 d(p) + m 2 d(p)) 2 — 4m 2 d(p) , (121) 

B(p, m, m) — l\ 



Ll(p,m) = -log 
LL(p,m 1 ,m 2 ) = log 



B(p, m, m) + 1 J ' 
1 — (m 2 + m 2 )d(p) + B(p, m 1 ,m 2 ) 
1 — (mf + m 2 )d{p) — B(p, mi, m 2 ) 



B.2.1 Result for the matrix elements of 4 in Full QCD with General kine- 
matics: 

We start with the matrix elements of O4 which we write in the form 

2 

(7r + 7r°|04|^ + > = Of c [1 + J^plzf] + 02 + 0\ + 0? d + counterterms , (122) 

where the counterterms are given in eq. (5). The two-pion state is symmetrised over the 
two-momenta p\ and p 2 to ensure that it has isospin 2. I z f is given in eq. (107) and we 
now present the remaining ingredients of the right-hand side of eq. (122). 

or = * {27 ' 1) x{3V2(PK-Pi+PK-P2 + 2p 1 - P2 )} ^ 

fnfK 
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a (27,i) i g l^fj (i2 Pl ■ p 2 + 8p K • Pi + %Pk ■ P2 ~ m 2 K ) 
a (27,i) bg (jpj [Am 2 K - ml) (2 Pl ■ p 2 + 3 (p K ■ Pl + p K ■ p 2 ) ) 
ml (86pi • p 2 + 29p K • Pl + 29p x • p 2 + 4ml) 



n log ^ 



16^2/ 2 /. 2 /x^ 2 



(124) 



^ = 



(27,1) 



4v/2/ 2 / 2 /*rr 2 

X 

rv(27,l) 



(div(pi +p 2 ) - log(d(pi + p 2 )m 2 )) 
6(pi • P2) 2 + 3(pi • p 2 ) (pk • Pi + Pk ■ P2 - 2m 2 ) 
+ ml (-p K ■ pi - Pk ■ P2 + 2m 2 

12V2P Pf K n* X { 36(Pl ' ^ + (Pl ' V2) ( 18PK ' Pl + 18PK ' P2 ~ 5 ° m 9 
+ ml (-9p K ■ pi - 9p K ■ P2 + 22m 2 ) 
+ 4d(pi + p 2 )m 2 

x 8(pi -p 2 ) 2 (2 + 3rf( Pl +p 2 )ml) 

+ m 2 ^3^ • pi + 3px • P2 + 4m 2 (-5 + <6d(pi + p 2 )m 2 ) 
+ (Pi • Pi) (?>Pk ■ Pi + 3p K ■ p 2 + Ami (-1 + 12d(pi + P 2 )m 2 ) ) J 
x |Ll(pi + p 2 ,m w )^2 - 8d(pi + p 2 )m 7r 2 

2 ( - 1 + 2d(pi + p 2 )m 2 ^) (p K • + p K ■ p 2 - 2m 2 + 8d(pi + p 2 )m^ 
+ 2(pi • p 2 ) 2 (3 + 6d(pi + p 2 )m 2 + 8d(pi + p 2 fm^j 
+ (Pi ■ P2) ( "ipK ■ P2 - 6m 2 + 2d( P i + p 2 )pk ■ P2ml + ?>2d(pi + p 2 ) 2 m^ 



(27,1) 



s/ 2 / 2 iw^ 



x 



+ 



(Pk-Pi) (3 + 2d( Pl + p 2 )ml) 



(125) 



and finally the monster 

n,(27,l) 



864 V2PPf K K 2 
x^3B(-p 1 + p K ,m v ,m K ) x LL{-p x + p K ,m v ,m K ) 
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(Pi • P2) [ - 9 (3p^ • pi + 8m^ - 2m; 

+ 2d(-p 1 + p K f (m 2 K - ml) {-2p K ■ p l + m\ + m 2 ) 

- + p K f (m\ - ml) 

x (^lm 4 K + AAm 2 K ml - 13m 4 + 18(p x • pi) (-5m^ + ml) 

- 3d(-pi + p K ) ( - 5m^ - \Am\ml + m 4 + 3(p K • Pi) (5m 2 K + ml 
+ (px • P2) [ - 9 (3pk • Pi + 10m^ - 4m^ 



- 2d(-pi + p K f{m 2 K - ml) 3 (-2p K ■p 1 + m 2 K + ml) 

+ 3d(-pi + p K ) (l9m 4 K - 2m 2 K m 2 n + m 4 + 9(p^ • pi) (-3m^ + ml) 

+ d(-pi + p K f (jn 2 K - ml) 

x (35rr4 + 44m^m 2 - 7m 4 + (p K ■ Pl ) (-78m|- + 6m 2 ) 



+2 log (d(-pi + p K )m K 

x 108d(-pi + p K f(m 2 K - ^) & (-2px ■ Pi + m 2 K + m 2 ^ 
- 2d(-p 1 + p^) 4 (2p^ -pi-m 2 K - ml) (m 2 K - m 2 ^ 

x (l21(p! • p 2 )m 2 K - 216m 4 K - 12l( Pl ■ p 2 )m 2 - 108m 2 K ml - 216m 4 . 

- V2\p K • p 2 (m 2 K - ml) + 270p^ • pi (m 2 K + m^) ^) 

+ 27^6(p^ • pi) 2 - 7(pi • p 2 )m 2 K + m 4 K + m 4 - (p K • p 2 ) {jm 2 K + 2ml) 

+ (Pk ■ Pi) (-Pi • P2 - 7p K • P2 + 3m^ + 3m 2 ) ^ 
+ 9d(-pi +px) 

x [ - ((Pi • P2)m A K ) + 21m^ - 16 (pi • p 2 )m\ml + 3m 4 K ml - (pi • p 2 )m 4 - 3m|- 

- 21m° + (p x • p 2 ) (-23m^ - 8m^m^ + 49m 4 ) 

+ 18(p x • Pi) (m 4 K - ml (2pi • p 2 - 2p^ • p 2 + m 2 ) ^ ] 
- 3d(-pi +p K ) 2 

x [ - 133(pi • p 2 )m 6 K + YJ\m\ + 138(pi • p 2 )m 4 K ml - lll(pi • p 2 )m^m 4 

- 110(pi -p 2 )m^ - 171m® - 90(p^ -pi) 2 (m 4 K - m 4 ) 

+ (Pk ■ Pi) (31m 6 K - \2m 4 K ml - 51m 2 K m 4 n + 248m^) 

+ 3(pty • Pi) ( - 51tt& + ll(pi • p 2 )m 4 n + 51m° - m 4 K (l9pi • p 2 + 15m 2 ) 
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+ (pk ■ p-i) (^3m% — 56m 2 K ml — 59m^j 

+ 5m 2 K (l6(pi • p 2 )ml + 3m*) ^ ] 

+ d(-ft+^) 3 ( m ^- m ^) 

x [ - 593(pi • p 2 )m x + 675m|- - 123(pi • p 2 )m^m^ + 324m^m^ 

+ 153(pi • p 2 )m^m 4 r + 162m 4 K m 4 n + 563(pi ■ p 2 )m\ + 324m^m^ + 675m* 
+ 540 (p^ • Pi ) 2 (m^ + m 2 ) 2 

+ (p* • P2) (461m^ + 2ttm 4 K m\ - 21m x m* - 695m^) 

- (®Pk ■ Pi) (252m 6 K + 151(pi • p 2 )m* + 252m° + m x (-161pi • p 2 + 108m 

+ 39 (px • P 2 ) (3m^ + 2m 2 K m\ - 5m^j 
+ 2m 2 K (hp 1 ■ p 2 m 2 n + 54m*) ^ ] 
-log (d(-pi + Pk)™ 2 ) 
x - (p K ■ P2) 

x [27 (3p x • Pl + lOm^ - Ami ) 

+ 2d(-pi + p K f (m\ - m 2 ) (-2p K • p l + m\ + m 2 ) 

- 9d(-pi + px) (25m 4 K + l6m 2 K ml - 5m* + 6(pk • pi) (Vm 2 

- d(-p 1 + PK) 3 (jn 2 K - m 2 ) (j7m 4 K + 80m^m 2 - 13m* 

- 18(px • pi) (9m 2 K - m 2 ) 
+ 3d(-pi + p K f (82m 6 K + l86m 4 K m 2 n - Wm 2 K m 4 + 8ml 

+ 3(p K -p 1 ) {-77m 4 K + 4m 2 K ml + m*) 

+ (Pi • P2) 

x [ - 27 (3p K ■ Pl + 8m 2 K - 2m 2 ) 

+ 2d(-pi + px) 4 (m^ - m 2 ) (-2p K • pi + m^ + m^j 

- d(-pi + Pk) 3 {tti 2 k - m 2 ) ^83m^ + 80m^m^ - 19m* 

- 6(pk ■ Pi) (29m x - 5m 2 ) 

- 9d(-pi + p^) ( Q7m 4 K - 32m 2 K m\ + m* + Q(p K ■ pi) (7m 2 K - m 2 ) 



4 



i) 
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+ 3d(-pi + p K f ( - 3{pk ■ Pi) (ll9m 4 K - 52m 2 K m 2 n + 5m 4 ) 

+ 2 (74rr4 + 8Am A K ml - ^7m 2 K m\ + 7m°) 
-27 \og(d(- Pl + p K )ml) 

12(p K ■ pxf + ll{p K ■ Pi)(p K ■ P2) ~ 6(pk ■ Pi)m 2 K + A(p K ■ p 2 )m 2 K - 2m% 



x 



3 



- Q(Pk ■ Pi)"^ + 8(p K ■ p 2 )m 7T - 2777^ 
+ 8d(-pi + p K f{m 2 K ~ "1^) (~2px • Pi + «4 + m 2 ) 
+ (Pi • P2) (~Px • Pi + 6m^ + 2m^) 

- 2d(-pi + p^) 4 x (2p K -pi-m 2 K - m 2 ) x (m 2 , - m 2 ) 

x [9(pi ■ p 2 )m 2 K - lQm 4 K - 9(pi • p 2 )m 2 - 8m 2 K m\ - 16m 4 . 

- 9(pk ■ P2) (m 2 K - m 2 ) + 20(p K ■ pi) (m 2 K + m 2 ) ] 
+ d(-pi +p K ) 

x [ - 23(j»i • P2)'m 4 K + 14777^ + ^ m K m l - (p 1 ■ p 2 )m 4 - 2m 2 K m 4 

- Urn* + (p K ■ P2) (-77T7 4 , + 31777 4 ) 

+ 2(p K -pi)( - 7(pi -p 2 )m 2 K + 6m A K - -p 2 )m 2 - 6m 4 . 

+ (PX"P 2 ) (7777^ + 117772)^) ] 

+ d(- Pl +p K f (m 2 K -ml) 

x [ - 47(pi • p 2 )m|- + 50tt& - 9(pi • p 2 )r?4m 2 + 24774m 2 + 15(pi • p 2 )m^m 4 

+ 127T7^777 4 + 4l(pi • p 2 ) 777^ + 24m^7T7^ + 50777^ + 40(j» X • Plf {m 2 K + 777 2 

+ (Pk ■ P2) (37777^ + lgm^m^ - 5m 2 K m 4 - 51m®) 

- 2(px • Pi) ^567?4 + 33pi • p 2 m 4 + 56m® + 37?4 (-13pi • p 2 + 8m 2 ) 

+ {p K ■ p 2 ) (29m 4 , + 14777^7772 _ 43m 4^ 
+ 6m 2 K ((Pi • P2)m\ + 4m 4 ) ^ ] 

- tZ(-pi +Pty) 2 

X [-20K.p 1 ) 2 (777 4 f - 777t) 

+ (Px • Pi) ( - 34m^ + 9pi • p 2 m 4 + 34777® + m 4 , (27pi • p 2 - 10m 2 ) 
+ 3(px-p 2 ) (< - 127772,7772 _ 13m 4) 

+ 2777^ ((I8pi • p 2 )7772 + 5m: 
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+ 2 x ( - 23(pi • p 2 )m K + 19m^ + 6(pi • p 2 )m 4 K ml - 6(pi • p 2 )m 2 K m 4 
- 13(pi • p 2 )m£ - 19m* 

+ (p* ' P2) (8m^ + 9m>^ - 9m^m^ + 28m°) ) ] 
+27B(-p 1 + p K ,m 7r ,m K ) x LL(-p x + p K ,m 7T ,m K ) 

x 12(pjr • pi) 2 - 11 (pk • Pi){pk ■ P2) + G(pk ■ Pi)m 2 K - A(p K ■ p 2 )m 2 K + 2m 4 K 
+ 6(Pk ■ Pi)™ 2 - 8(p K ■ p 2 )m\ + 2m% 
+ 8d(-pi + pK) 4 (m 2 K - m 2 ) (-2p K ■ p 1 + m 2 K + m 2 ) 

- 6d(-pi + p^) 3 (2p K -px-m 2 K - m 2 ) (m 2 K - m 2 ) 
x ^3(pi • p^m^ - 4m^ - 3(pi • p 2 )m 2 - Ami 

- 3(Pk ■ P2) (m 2 K - m l) + ^{Vk ■ Pi) (m\ + m 2 ) ^ 

+ (Pi • P2) (pk ■ Pi - 2 (3m^ + m 2 ) ^ 
+ + p^) 2 

x [ - 29(pi • p 2 )m K + 26m 8 K + 9(pi • p 2 )m 4 K m\ - 8m K m 2 K - 3(pi • p 2 )m 2 K m 4 

— Am 4 K m 4 + 23(pi • p 2 )m\ — 8m 2 K m\ + 26m® 

- 8(p^ • pi) 2 (m^ - Qm 2 K m 2 n + m 4 } 

+ (Pk ■ P2) (l9m% + m K ml + 13m 2 K m 4 - 33m^) 



4 

m K 



-2(p K -p 1 )(20m% + 5mi (?,p 1 ■ p 2 + Am n 

+ (Pk ■ P2) (llm^ + \Am\m\ - 25m*) 
+ m^ (6(pi • p 2 )m 2 - 4m^) ^ ] 

+ tZ(-pi + p*r) 



(21pi • p 2 + 4m 2 



x [17(pi • p 2 )m 4 K - 12m K - 8(p l ■ p 2 )m\m\ - 3(pi • p 2 )m\ - 12m^ 



+ \2{p K ■ pi) 2 (m 2 K + m 2 ) + (pk • p 2 ) (3m^ - \2m\m\ + 23m^ 
+ (Px • Pi) ((pk • p 2 ) (-25m^ + 11m 2 ) 

- 3 (2m^ + 7(pi • p 2 )m 2 + 2m^ - m 2 ^ (5pi • p 2 + 4m 2 )) 



+ 



(Pi ^P 2 )| 



a 



(27,1) 



X 



72d(-pi + PK) A {rn 2 K - m 2 ) (-2p^ • pi + m 2 ^ + m 2 ) 



144 V2f 2 f 2 f K 7r 2 

- 72d(-p 2 + p K ) 4 (m 2 K - m 2 ) (-2p K • p 2 + + m 2 ) ' 
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+ Ad(-pi + p K f fax -Pi-m K - m 2 ) (m 2 K - m 2 .) 

x ^41(pi • p2)m K - 63m K - 41(p 1 • p 2 )m\ - \8m\ml 

- 63ml - 41(p K ■ pa) (m^ - ml) + 72(p K • p x ) (m\ + m*) ^ 
+ 4d(-p 2 + p K f {%Pk -P2-m K - m 2 w ) (m 2 K - m 2 ^ 

x ^41(p! • P2)m 2 K - 63m^ - 41(pi • p 2 )m^ - \8m\m\ 

- 63ml - 41(p K ■ p x ) (m 2 K - ml) + 72{p K ■ p 2 ) (m 2 K + ml) ^ 
- 18div(-p 2 + Pk) (g(pk • P2) 2 - 7pi • P2m 2 K + m K + ml 

-Pk-Pi (jPK ■ P2 + 7m 2 K + 2ml) + p K ■ p 2 {-pi ■ p 2 + 3m 2 K + 3ml) ^ 

- 3 x [ 36 (2 + div(-pi + p K ) ) (p K • Pi) 2 + 72{p K ■ P2) 2 - 278( Pl ■ p 2 )m 2 K 
- 42div(-p! + p K ){p\ • P2)m 2 K + 42m K + 6div(-pi + p K )m K 
+ (lppi • p 2 )ml + Am 2 K ml + 42ml + 6div(-pi + p K )ml 

+ (Pk ■ P2) ( - 28 Pl • p 2 - (71 + 42div(- Pl + p K )) m\ 
+ (l-12div(-p 1 +p^))m^ 

+ (Pk -Pi)(^- 28pi • p 2 - 6div(-pi + p K )(pi • P2) 
- 2 (100 + 21div(-pi + p K )) (pk ■ P2) 

+ (-71 + 18div(-p! + p K )) m 2 K + ml + 18div(- Pl + p K )m^j ] 

+ 3d(-p 2 +p K ) 

x [ - 89(pi • p 2 )m K + 65m 6 K - 24{ Pl ■ p 2 )m 2 K ml + 15m K m 2 n + 41{ Pl ■ p 2 )m A n 
+ \hm\ml + 65ml - 64(p K • p 2 ) 2 (m 2 K + ml) 

-2(Pk- P2) (7m K - 16pi • p 2 m\ + 7m 4 n + 2m 2 K (-10pi • p 2 + m 2 ) ) 
+ (Pk ■ Pi) (-31m K + 24m 2 K ml - 113m* + 8p K • P2 (l3m 2 K + 2ml) ) ] 
+ 3d(-pi +p K ) 

x [ — 89(p 1 • p 2 )m K + 65m K — 24(p x • p 2 )m 2 K ml + Yhm K m\ 

+ 41(pi • p 2 )ml + \hm\ml + 6hm% - 64{p K • Pi) 2 (m 2 K + m 2 ^ 
+ (Pk ■ P2) (~31m K + 24m 2 K ml - 113m*) 

+ [2p K • Pi) - 7m K + 16pi • p 2 m\ - 7ml - 2m 2 K (-10pi • p 2 + m 2 ) 
+ {Pk ■ P2) (52m 2 K + 8ml) ) ] 
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- 2d{~p 2 +p K ) 2 

x [ — 202(px • p 2 )m 6 K + 165ml- + 66(pi • p 2 )m 2 K m 4 — 42m 4 K m\ 
+ 136(pi • p 2 )m\ + lQ5ml + 48(px • P2) 2 (jn% + A.m\m\ + m 4 

+2{p K ■ pi) [m 2 K - ml) (j7m 4 K + \^\m\m\ + 92m^ - 30p K ■ p 2 {^rn 2 K + 5ml) 

- 6(p^ ■ p 2 ) ^53m|- + 34(pi • p 2 )m^. + 53m^ - m 4 K (56pi ■ p 2 + 5ml) 

+ m 2 K (22 Pl -p 2 ml - 5mf) ) ] 

-2d(- Pl +p K ) 2 

x [ - 202( Pl • p 2 )m^ + 165m^ + 66(pi • p 2 )m 2 K m 4 - 42m 4 K m 4 n 
+ 136( Pl • p 2 )m6 + 165m« + 48(p x • Pl ) 2 (m 4 K + 4m 2 K ml + m 4 ) 

+ 2(p K ■ p 2 ) (77m% + 2Am 4 K ml - 9m 2 K ml - 92m^) 

- 6(p K ■ pi) ^53m^ + 34pi • p 2 m^ + 53m^ - m 4 K (56pi • p 2 + 5ml) 

+ m 2 K ((22 Pl -p 2 )ml-5m 4 n ) 
+ 10(p K ■ p 2 ) (^m 4 K + m 2 K ml - 5m 4 w ) J } 



B.2.2 Result for the matrix elements of (9 7 8 in Full QCD with General Kine- 
matics: 

For the EWP operators the corresponding expressions are: 

2 

(7r + 7r°|0 7i8 |^+> = 0\J [1 + j^phf] + 7 % + O b 7>8 + Off + counterterms , (127) 
where the counterterms are given in eq. (122), I z f is given in eq. (107), 

(128) 



37 <s ' 8) log j m\ 77M log (f) ml 7 <8 ' 8 > log [jf ) ( + < 

^7,8 , it: r n „ 9 „ _ 9 „ p: r0 „ 9 „ _ 9 



4V2f 2 f 2 f K ^ sV2f 2 f 2 f K ^ 24V2f 2 f 2 f K ^ 

(129) 

Qb 7 (s,8) (div(p! + p 2 ) - log(^(Pi + P2)ml)) (3pi • p 2 - rrffi 

7,8 ^ " QV2f 2 f 2 f K 7T 2 

7 (8,8) [ m 2 (_ 3 + 4rf ( pi + p2 ) m 2) + ( pi . p2 ) (g + 4^( pi + P2 )ml)} 

6V2f 2 f 2 f K rr 2 
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(126) 



and 



O c + d 



x 



7 ( 8 ' 8 )Ll(p! + p 2 , m«)j2 - 8d(pi +p 2 )ml 



x 



m 



2 (-1 + 2d(px + p 2 )m 2 ) + ( Pl ■ pa) (3 + 2d(p! 



(130) 



7 



(8,8) 



864 V2/ 2 /2fr7r 2 

I - 3 + p x , m,, m K ) x LL{-p x + p*r, m,,, %) 

x — 27p^ • pi — 6m|- + 6 m 2 . 
+ d(-pi + p x ) 2 (m 2 ? - m 2 ) (-2 p K ■ p 1 + m\ + m 2 ) 
+ <i(-pi+px) (-7mjf 4 - 16 m 2 K ml + 5 m 4 + 3p K • pi (Vm^-m 2 )) 
+3 5(-p 2 + Pk, m,,, m K ) x LL(-p 2 + p*r, m,,, m^) 
x - 27p^ • p 2 - 6mjf + 6 m 2 

+ d(-p 2 + p K f (rn 2 K - m 2 ) (-2 Pifp 2 + m^ + "4) 
+ <i(-p 2 +px) (-7m^ 4 - 16 m 2 K ml + 5 m 4 . + 3 p K ■ Vi {j ' m K ~ m l)) 
+ log (d(-p! + p K )mfj 
x 9 (9p^-pi + 2m^-2m 2 ) 

- 9 d(-pi + p x ) (m^ - m 2 ) (-2 p x • pi + 7 m^ - m 2 ) 
+ d(-p 1 +p K f (rn 2 K -m 2 ) 3 (-2p x ■ p 1 + m 2 K + m 2 ) 

- + px) 2 (jn 2 K - m 2 ) (28 m 4 ^ + 34 m 2 K m 2 - 8mJ + 9 px • Pi (-7 m 2 ^ + m 2 
+ log (d(-p 2 + p K )mty 

x 9 (9p^ -P2 + 2 771I- -2m 2 ) 

- 9 d(-p 2 + p K ) (rn 2 K - m 2 ) (-2 p K ■ p 2 + 7 m\ - m 2 ) 
+ rf(-p 2 +px) 3 (m^ -m 2 ) 3 (-2p x -p 2 + m^ + m 2 ) 

- d(-p 2 + Px) 2 - ^ 2 ) (28 m A K + 34 m^ m 2 - 8 m 4 + 9 p K ■ p 2 (-7 m|- + m 2 
+9 5(-pi + p x , m v , m K ) x LL(-p x + p^, m,,, mjf) 

21p K -p 1 -4 (m^ + m 2 ) 

+ 7 + p K f [m 2 K - m l) 2 (-2 Pk ■ Pi + m 2 K + m 2 ) 
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x 



+ d(-pi + p K ) (-3 m K - 8 m 2 K m\ - 3mJ + 7p K ■ pi (m 2 K 
+9 B{-p 2 + p K , m v , m K ) x LL(-p 2 + p K , m v , m K ) 
x 21 pa: ■ p 2 -4 (m^+mj) 

+ 7 d(-p 2 +Pk) 2 (jn 2 K - ml) (-2 p K ■ p 2 + m 2 K + m 2 ^ 



+ c?(-P2+Px) (~3m% - 8m 2 K ml - 3ml + 7 p K ■ p 2 (m 2 K + ml)) 
+6 x - 22 d(-pi + p K f (m 2 K - ml) (-2 p K ■ p 1 + m 2 K + ml) 

- 22 d(-p 2 + p K f (m 2 K -m 2 ) (~2p K ■ p 2 +m 2 K + ml) 
+ 6 [-15 (2 + div(-pi + p K )) p K ■ Pi ~ 15 (2 + div(-p 2 + p K )) p K ■ p 2 

- 3 m 2 K + div(-pi + p K ) m 2 K + div(-p 2 + pk) m 2 K + 9 ml 
+ 3 div(-pi + p K ) m 2 w + 3 div(-p 2 + p K ) ml] 
+ d(-p 1 +p K ) (37 m 4 K + 70 m 2 K ml + 13 m 4 n - 12 p K - Pl (jm 2 K + 3ml)) 

+ d(-p 2 +p K ) (37 m K + 70 m 2 K ml + 13 ml - 12 p K ■ p 2 (7mj s - + 3mjj) 

+2 log (d(-pi + PK)m 2 K )) 

x 9 d(-pi + pat) (m^ - m^) (6 p x • Pi + 3 m|- - m£) 

+ 31 d(- Pl + pk? (m 2 K - mlf (-2 p K -pi+m 2 K + ml) 
-9 (-15pA--pi+m^ + 3m^) 

-d(- Pl +p K ) 2 (m 2 K -ml) (31 m 4 K + 82 m 2 K ml + 49 m%- 9 p K • Pi (7m|- + 11 to')) 

+2 log ((d(-p 2 + PK)m 2 K ) 

x 9 d(-p 2 + p K ) {rn 2 K - m 2 w ) (q p K ■ p 2 + 3 m 2 K - m 2 n ) 

+ 31 d(-p 2 + pk) 3 (m 2 K - ml) 3 (-2 p K ■ p 2 + m\ + m 2 w ) 

- 9 (-15 p K ■ P2 + m\ + 3 ml) 

-d{-p 2 +p K f (m 2 K -ml) (31 m 4 K + 82 m 2 K ml + 49 m%- 9 p K ■ p 2 (jm 2 K + 11 ml)) 
-9 log (d(-pi + p K )ml) 
x - 21p K -p! + 4 (to|- + m 2 n ) 

- d(-pi + pat) (m 2 K - ml) (-14 p^ • pi + to|- + to') 
+ 7 rf(-pi + p K ) 3 (m 2 K - ml) 3 (-2 p K ■ pi + m 2 K + ml) 
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+d(-pi + p K f (21 Pk ■ Pi (m% - m^j - 2 (5 m\ + 6 m A K ml - 6 m 2 K m\ - 5 m®)) 

-9 log (d(-z>2 +P*)"4) 
x -21p K -p 2 + A (m 2 K + ml) 
- d(-p 2 + px) (m^ - ml) {-I^Pk ■P2 + m 2 K + ml) 
+ 7 <Z(-p2 + P^) 3 (m^ - m 2 n ) 3 (-2 p K • p 2 + m\ + m 2 ) 
+d(-p 2 + Pk) 2 (21px • ("^ - ra*) - 2 (5 + 6 m 4 K ml - 6 m 4 . - 5 ^)) | • 



:i3i^ 



C Results in Quenched QCD 

Finally we present the results for matrix elements in quenched QCD. The results for general 
kinematics, i.e. for arbitrary quark masses and momenta are too lengthy to be exhibited 
in this paper; we present them on the web site [23]. In the following subsection we present 
the results for the matrix elements with the SPQR kinematics. 

In the quenched approximation, the integral 7 , which corresponds to the one-loop contri- 
butions from the renormalization of the mesonic wave-functions and from the replacement 
of the factor 1/f 3 which appears at lowest-order in the chiral expansion, by is 
given by (y = ml/m 2 K ): 



2ml 
9m 2 K 



log 



2(1 -y) °\2-y 



2a 



i + !^4iog| 



2(1 -y) °V 2 -», 



(132) 



As for full QCD, the results in this appendix are defined in terms of the variables, y = 
ml/m 2 Kl z = m n /m K and u> = E n /m w , where E n is the energy of the pion whose momentum 
is (in general) not equal to zero. 

C.l K — > 7T + 7r° Decays with SPQR Kinematics in Quenched QCD 

C.l.l Result for the matrix elements of O4 in Quenched QCD with SPQR 
kinematics 



The matrix elements of O4 with the SPQR kinematics in the quenched approximation is 
given by 



(7T + 7r°|e>4|i^ 



-6^77^(27,1) f z(l+u + 2z) f m\ 



f K f 2 



2u 



16tt 2 / : 



in 



K 



16tt 2 / 2 



(i q a + n + it 



c+d f' 



(133) 



7 We introduce a superscript q in this appendix to denote the quenched approximation. 
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where 1^ is given in eq. (132) and the I^c+d are as follows: 
11 = l^{(-2z 2 +u-2z(l + u;))\og(z 2 ) 



and 



1 c+d 



+ (2z 2 -uj + 2z(l + uj) + 2z 3 (l + uj) + z 4 (Q + u)) log (^pj j 

a (3 + 2z + 3uj) (2 - 2z 2 + z 2 (-2 + z 2 )\og (-1 + J,)) 

18^ 3 (-1 + ^ 2 )cj 
m 2 (3 + 2z + 3cj) (-2 + 2z 2 + log (-1 + J,)) 
18z(— 1 + 2 2 )u;m 2 r 



4 9 



3w 2 



+ 



){ 



3(1 + 2z + u) + 3W^— ^(1 + 2^ + w)log 



1+U1 



1+UJ 



1+u; 



-3 - 2cj + cj 2 - 2^(3 - 3c; + cu 2 ) log 



m 



/i 



-1 



x 



-3 



18(— 1 + z)uj 2 {-2z + uj + z 2 ujf / 
IOzuj 4 -2uj 5 + 2z 9 uj 2 (3 + cj 2 ) -2^ 2 cu 3 (9-6uj + 5uj 2 ) 

cu 2 (-22 + 56 uj - 41 cj 2 + uj 3 ) - * 8 cu 2 (l2 + 42 - 9 uj 2 + cu 3 ) 
+ (-32 + 90cj-64cj 2 + 15cj 3 -17c/) 
+ z 7 (-8 + 24cu + 56cu 2 - 18 cj 3 + 55 cj 4 - 3 cj 5 ) 
+ 2/ (12-18cj + 28cj 2 -33cj 3 + 36cj 4 + 6cj 5 ) 

- 2 * 6 (8 + 14 uj - 23 cj 2 + 48 uj 3 + 13 uj 5 ) log(z 2 ) 
6 (-I + 2) * 3 Vl -uj 2 [3 2 5 cu 2 + 3 cj 2 (-1 + 2uj) ~3z a uj 2 (l + 2u) 



-3zuj (-4 + 6CJ + CJ 3 ) +Qz 2 (-2 + 3uj-uj 2 + uj 3 ) 
+ z 3 (-4 + 12cj-7cc> 2 + 5c</) 



log 



/ 1 - VT^J 2 " 



(134) 



(135) 



+ (-l + z) (-2 z + uj + z 2 uj) 2z 2 ^-§uj 2 (l -4cj + cj 2 ) + z 5 uj 2 (7 + 9uj + 2uj 
+ zuj (36 - 128cj + 45cj 2 - 7cu 3 ) - z 4 uj (-2 + 45 uj + 2 cj 2 + 9 uj 3 ) 
-2z 2 (18-64cj + 45cj 2 -26cj 3 + 9cj 4 ) + z 3 (-8 + 72cj- IOIcj 2 + 54cj 3 + 19c/) ) 
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+3 (-2 



Z + UO + z uo 



2uo 2 + 2 z 4 uo 2 + 2 zuo (1 + uo) + z 3 (2 + 3uo + uo 2 ) 



+2z 2 (3-5cj + 3cu 2 )")log 



a 



18(-1 + zf(l + zfuo\-2z + uo + z 2 uf i 
1 



X 



2zuo - 2z^uo -uo 2 + z^uo 2 - z 2 (-l + cu 2 ) ) V( 1 + ^) 3 (1 + ^) x 

— ^(-1+^)3(1+^) 



2zu; - 2z 3 uo -uo 2 + z 4 uo 2 - z 2 (-l + uo 2 ) 



X 



z[-12 uj 4 + 4z 11 uj 4 + z w uo 3 (-6 + 9uo - 5 uo 2 ) - 2 z uo 3 (-23 + 8uo + uo 2 ) 

+ 2z 2 uo 2 (-45 + 61cj + 9w 2 + cj 3 ) + z 9 cj 2 (l8 - 72uo - 3uo 2 + 5cj 3 ) 
-2z 3 cj (-27+ 123cj-60cj 2 + 54cj 3 + 14cu 4 ) 
f42 -92w + 15cj 2 -16cj 3 + 17cu z 



+ 2z 8 uo (-10 + 69cj-45cj 2 + 18cj 3 + 17cj 4 ) 

-* 6 u; (78 + 100cj + 42cu 2 + 39cj 3 + 63cj 4 ) 

+ 2z 4 (-2 + 85 - 108 cj 2 + 126 uo 3 - 54 c/ + 14 cj 5 ) 

+ 2 5 (4 + 126 uo - 192 cj 2 + 204 uo 3 + 19 uo 4 + 63 cj 5 )] log (-1 + — ) 

-2(l + z)o;^(-l + z) x 

[- (uo 3 (1+uo)) + 3z 9 uo 2 (-4-3cj + cj 2 ) + zuo 2 (6 + 27 uo + uo 2 ) 
+ z s uo (24 + 66 uo + 25 uo 2 - 3 uo 3 ) - z 2 uo (l2 + 107cj + 17cj 2 + 22 uo 3 ) 
-z 1 'uo (132 + 79cj + 39cj 2 + 22cj 3 ) -Az b (3 - 8uo - 8uo 2 + 14cj 3 + uo 4 ) 
+ 2z 4 (-46 - 8uo - 64uo 2 - 35uo 3 + 2t 
+ z 6 (72 + 16cj + 109cj 2 + 67cu 3 + 22( 

+ z 3 (8 + 160cj + 41c; 2 + 97cj 3 + 22cu 4 ) 

+ z(l + z) 6 uo 3 + 3 cj 2 (-3 - 2 + uo 2 ) +zuo 2 (-23 - 3 uo + 2 cu 2 ) 

+ 3 z 7 uo (6 + 21 uo + 4uo 2 - 3 uo 3 ) + z 6 uo (-126 - 89 - 33 uo 2 + 2 uo 3 ) 
-3z 2 uo (-\A + 7uo + 3uo 2 + 2uo 3 ) + 3 z 5 (l6 + 40 + 47 cj 2 + 19 uo 3 + 6 uo 4 ) 
+ z 3 (-24 + 54cu + 11cj 2 + 78cj 3 + 13cj 4 ) 



:oo 



:to 
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- z A (48 + 36 uj + 145 uj 2 + 78 uj 3 + 23 uj 4 ) log(z 2 

+ 3(1 + ^) 2 (-4 + 10 2: -3z 2 -7 z 3 + Az 4 )uj (-2z + uj + z 2 uj Y x 
2^cu-2z 3 cj-cu 2 + 2 4 cj 2 -2 2 (-1+w 2 ) x 

'l + z-z 2 - ^J(-l + z) 3 (1 + z) 



log 



l + z-z 2 + ^{-l + zf (1 + z). 



- (-1 + z)z 2 J(-l + z) 3 (l + z) X 



- 6cj 4 - Qz w uj 4 + 27 z 9 uj 5 + 9uj 6 -9zuj 3 (-2 + 3cj 2 ) 



z 6 oj 2 



-38 + 67 uj 2 + 55 uj 4 ) 



-3z 7 u (a~3uj 2 + 7uj 4 ) 
-3z 8 uj 2 (-6 + 10cj 2 + 7 to 4 ) + 
+ z 3 c; (2 -67cj 2 + 113 c/) + z 5 (lO cj + 40 u 3 - 92 a/) 
+ ^ 2 (-10cj 2 + 29cj 4 -25cc> 6 ) 

/l - y /T=u s \ 



log 



+ z 4 (4 + 18 cj 2 - 42 c/ -22 c/) 

+ ^^2zuj - 2z 3 cj - c/ + /c/ - z 2 (-\ + uj 2 ) x 

- 12c/ + 6z n i/ (-3 + c/) + 22zw 4 (2 + uj 2 ) 



Kl + ^T^uT 1 



r 2 , , 3 
Z UJ 



(38 + 73 c/) 

-/cj (-12 + 107 c/ + c/) -2z 3 uj 2 (-8 + 15c/ + 5c/ ) 
+ 2 / c/ (48 + 55 to 2 + 35 c/) 

-2/cj (7 + 50 c/) + z 4 uj (2 + 97c/ + 177 c/) + z 10 (48 c/ + 9 

+ z 9 (-42 c/ + 11c/ -35 c/ ) -z 5 (4 + 58c/ + 129c/ + 49 

( z + uj - z 2 uj - J2zuj - 2z 3 uj -uj 2 + z a uj 2 - z 2 (-l + uj 2 ) ' 

log V , 

\z + uj-z 2 uj+ ^2zuj - 2z 3 uj -uj 2 + z a uj 2 - z 2 {-l + uj 2 



UJ 



+ m x 



18(-l + zf(l + z) 2 uj\-2z + uj + z 2 ujfm\ 
1 



x 



2zuj - 2z 3 uj -uj 2 + z a uj 2 - z 2 (- 
-1 

z 



— ^(-l + z) 3 (l+z) 



22cj - 2/cu - uj 2 + z 4 uj 2 - z 2 (-l + uj 2 ) 



x 



[- (uj 4 (1 + uj)) -3z 9 uj 3 (-3-2cj + c/) +zuj 3 (6 + 15cu + c/) 
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+z 8 uj 2 (-12-45cu + 4cu 2 + 3cj 3 ) - z 2 uj 2 (l2 + 49 uj + 2 uj 2 + 5 cj 3 ) 
+ z 3 cu (8 + 74cj-35cj 2 + 18cj 3 + 5cj 4 ) 
+ 2 4 uj (-42 + 106 - 81 uj 2 + 33 uj 3 + 12 cj 4 ) 
+ /cj (6 + 84cj-75cj 2 + 30cj 3 + 13cj 4 ) 
+ (4 - 86 uj + 110 cj 2 - 141 uj 3 + 27 c/ - 12 cj 5 ) 

-2 6 (4 + 54cj-154cj 2 + 93cj 3 -38cj 4 + 13cu 5 )] x log (-1 + — ) 



+ 2(1 + 2> 



(2 + 5cj + cj 2 )) 



(-1 + z) (-2 2 + w + z 2 uj x 
(4 z A u - uj (1 + u) - z 3 (3 + uj 2 )+z 2 (-5-2uj + u 
-z (1 + z) ( - 6 cj 3 - 3 ^ 5 cj 2 (-10 - 3 uj + cj 2 ) + 2 6 cj 2 (-4 - 3 uj + cj 2 ) 
+2V (20 + 9cj + 7cj 2 ) -2z 4 u (l8 + 7cj + 18 uj 2 + 2 cj 3 ) 
-3z 2 cj (12 + 2c; + 9c«; 2 + 7cj 3 ) + ^ 3 (24 + 46 uj 2 + 30 cu 3 + 20 uj 4 ) ) log( 



+ 3(1 + z) 2 (l - 3^ + 2^ 2 ) uj (-2 z + uj + z 2 uj) x 
(2 2cj-2 2 3 cj-cj 2 + 2 4 cj 2 -^ 2 (-l+cu 2 ))log 
- (-l + z)^(-l + z) 3 (l + z) x 



l + 2-z 2 - v/(-l + ^) (1 + *) 



l + 2-2 2 + + ^) (1+*) 



v / T^2[6,2 10 c/ + 9^ 5 -3a/ + 3,2 2 cc; 4 (-3 + 5uj 2 ) - 3 z 9 w 3 (2 + 7cu 2 ) 

+2^ 4 c; 2 (-1 + 17cj 2 + 5cj 4 ) + z 8 uj 2 (-10 + 32cj 2 + 11 uj 4 ) 

+z 7 uj (l4- 13 cu 2 + 23 cj 4 ) +2z 5 uj (-7 + 2uj 2 + 26c/) + (l5cu 3 - 63cj 5 ) 

-2 6 (4-24cu 2 + 75cj 4 + 29cj 



log 



l-VT 
1 + VT 



+ -v/2^ - 2^ 3 cj - cj 2 + 2 4 u; 2 



-1 +cu 2 ) x 



[3cj 5 + 2z 9 w 4 -4 + u> 2 -2w 4 4 + 5cj 2 +/cj 3 14 + 13^ 



7 4 

2; a; 



(23 + 13 cu 2 ) +z 2 uj 3 (-1 + 37cj 2 ) -z 3 uj 2 (4 + 3 uj 2 + 5 uj 4 ) 
-2z 4 uj (-5 + 8cj 2 + 24c/) + z 6 (-10 uj + 3 uj 3 - 5 cj 5 ) 
+ 2 5 (4-8cj 2 + 54c/ + 22cj 6 )] 



xlog 



z + — 2; a; 



2zuj - 2z 3 uj -uj 2 + z 4 u 2 



-1 +UJ< 



z + uj-z 2 uj + J2zuj - 2z 3 uj -uj 2 + z 4 uj 2 - z 2 (-l + uj 2 



(136) 
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C.1.2 Results for the matrix elements of (9 7 8 in quenched QCD with SPQR 
kinematics 



The matrix elements of the EWP operators are given by: 



<7rV>|0 7>8 |*+> = |4^ (8 ' 8) 



1 + 



m 



K 



lQir 2 f 



(jzf + Ja + Jb + Jt+d 



(137) 



where I z f, the one-loop contribution from the renormalization of the mesons' wavefunctions 
and the replacement of / 3 by fxfl is given in eq. (132) and the J a ,b,c+d are as follows: 



T q = 



(\og(z 2 ) + (-1 + z 2 )\og(^)) a (2 - 2z 2 + z 2 (-2 + r>g(-l + 



9(-l + z 2 ) 



2 z 2 (-2 + 2^ + log(-l + ^)) 
9(-l +z 2 )ml 

-2z 2 U + Sy^log (^jk ] + (-3 + c)log(5 



(138) 



3a; 



(139) 



and 



Jc+d — 



x 



3 (-1 + 2) o> (-2 z + o> + ^ 2 to) 

2 + zu (l + 2cu) + z 2 (-3-7cj + cj 2 ) + r 3 (3 + 2cj 2 )] log (V) 



or 



+3 (-I + 2) ^Vl-u; 2 lo£ 



fi-VT 



Vi + Vi-w 2 / 

+ (-1 + 2) (-2 ^ + + ^ 2 -3 2 (l+u) + (w + z 2 w + 3z + logfej J 



+ a 



\18(-1 + zf(l + ^) 2 cj(-2^ + uj + ^ 2 cj) 2 

^ZCJ - 2Z 3 LU -LU 2 + Z^LU 2 - Z 2 (-l + UJ 2 ) 



X 



-1 + zf(l + z) X 



2m> - 2z 3 u; - cj 2 + z 4 ^ 2 - z 2 (-l + cj) 2 



X 



z[z 6 (5 -7a;) + / (-5 + w) w 2 + 2 2 8 cj 3 - 3 a; 2 (1 + uj) + zuj 2 (15 + u) 
+ z 5 to (7 - 10 uj + 3 a; 2 ) + 2 (4 + 5 uj + 3 a> 2 ) 

(9 
-1 + — 

+ (-l + z 2 ) uj[- (z (-5 + u) u) -2uj 2 -Az 5 uj 2 + z 4 uj (13 + a;) 
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- (2 + z - 



-z 2 (2 - 7 uj + uj 2 ) -z 3 (lO + uj + 5 uj 2 ) 

+ 3z (1 + z) (uj (1+uj) + z 4 uj (1 + uj)-2zuj (3 + u) 

2 z 3 uj (3 + u) + 2 z 2 (2 + uj + 3 cj 2 ) ) log(^ 2 ) ] 

+ 2^ 4 ) (-2z + + z 2 u;) 2 x 
(-1+cj 2 )" x 



4^ 



2 ^3 



o o 3 2,4 2 2 

2zu; — 2z UJ — UJ +Z UJ —z 



l + z-z' 



v l + 2-2 2 + v /(-l + ^) 3 (1 + 2). 



log 

- 3(-l + ^) v /(-l + 2) 3 (l + ^) x 
./i ^[-2z 7 cu 3 + z 2 (2-5cu 2 ) +cu 2 (-2 + 3cu 2 ) + z 6 u> 2 (4 + 3cu 2 ) + 2z 3 u; (-2 + 7 

+z (4cu-6cu 3 ) -2^ 5 (cu + 2cu 3 ) (u; 2 -6u; 4 )]log ( |~ ^~^ ) 
+ ^2zuj-2z 3 uj-uj 2 + z 4 uj 2 -z 2 (-1+UJ 2 ) 



u 



x [22 7 cu 3 -2cu 4 + 2 5 w (2 + 3cj 2 ) -z 6 uj 2 (4 + 3cj 2 ) 
+ zuj (-2 + 7cj 2 ) +z 3 (2uj- Uuj 3 ) + z 2 (-2 + 5 cu 2 - 2 c/) 



(4 + 



(cu 2 + 7c/)] x 



log 



' z + uj — z 2 uj — ^2 zuj — 2 z 3 uj — uj 2 + z a uj 2 — z 2 (— 1 + uj' 
+ uj-z 2 uj + ^2zuj-2z 3 uj-uj 2 + z 4 uj 2 -z 2 (-1 + uj 2 ) # 



\18(-1 + ^) 5 (1 + z) 2 u;(-2z + + z 2 w) 2 m2 / X 

(2^-2^-, 2+ U2 + w 2)) *V(" 1 + *) 3 ( 1 + *)X 

^^(-l + ,) 3 (1 + z) 2zc; - 2, 3 cu - cu 2 + 2V - , 2 (-l + uj) 2 x 

2(-l + z 2 V[z(-5 + uj)uj + 2cj 2 + Az b uj 2 - A;(13 + w) 

+ z 2 {2 -7uj + uj 2 ) + z 3 (10 + uj + 5uj 2 )] 
+ [z (-5 + uj) uj 2 + 2uj 3 + 3 z 6 uj 2 (1 + uj) + z 5 uj (-6 - 15 uj + uj 2 ) 

+ z 2 uj (8-7uj + uj 2 ) +2z 3 (-3 + 2uj -7uj 2 + 4uj 3 ) 
+ 2z 4 (3 + 9^-5^ + 4^)] log + — ' 
- 3(-l + z 2 )u(-2 z + uj + z 2 ujf x 
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2 zoo - 2 z 3 uo - uo 2 + z 4 uj 2 - z 2 (-1+cj 2 ) 
log 



X 



1 + z-z 2 -^(-1 + zf (1 + z) 



l+z-z 2 + yj(-l + zf (1 + z), 
- 3(-l + Z )y/(-l+zf(l + z) X 

Vl -uj 2 [6 z 5 uj + 3 z 2 uj 2 + 2 z uj 3 + 2 ^ 7 uj 3 - uj 4 - z 6 uj 2 (6 + cj 2 ) 
-2^ 3 cu (2 + 3u 2 ) + z 4 (-2 + 5 u 2 + 2 cu 4 ) ] log ^i^L^j 

+ y/2zu-2z 3 u -uj 2 + z a uj 2 - z 2 (-1+cu 2 ) X 
[2 2 2 cj-32cu 2 + c; 3 + z 4 cu (-4 + w 2 ) 

- z 5 uj 2 (-2 + cj 2 ) + z 3 (2-uj 2 + uj 4 )} x 



z + uj-z 2 uj- J2zuj -2z 3 uj -uj 2 + z 4 uj 2 - z 2 (-1 + 



log V • (140) 

V^ + cj-^ 2 cj + ^2zuj-2z^uj-uj 2 + z a uj 2 -z 2 (-1 + uj 2 ) ) J J 
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